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Preface 


This  research  was  performed  to  investigate  site  parameters  affecting  natural  attenuation 
of  JP-4  in  groundwater  in  the  saturated  rone  and  to  suggest  general  site  conditions  which 
favor  natural  attenuation.  If  this  evaluation  is  of  some  assistance  from  either  a 
methodological  or  results  standpoint,  my  efforts  will  have  been  worthwhile. 

I  was  attracted  to  the  study  of  natural  attenuation  because  I  have  an  affinity  toward  the 
natural  -  -  natural  foods,  natural  energy  sources,  natural  healing,  etc.  This  remediation 
alternative  provided  a  natural,  common  sense  approach  to  a  serious  problem.  Although  it  may 
take  more  time  than  technologically  based  alternatives,  natural  attenuation  promises  to  be  a 
reasonable  solution  to  the  escalating  costs  of  ate  remediation.  In  becoming  more  familiar 
with  the  process,  I  found  that  many  authors  have  recognized  the  advantages  of  this  option  and 
ongoing  research  is  being  done  to  better  understand  the  process. 

I  wish  to  express  my  gratitude  to  Lt  Col  Michael  Shelley,  my  thesis  advisor,  and  Lt 
Col  Mark  Goltz,  member  of  the  thesis  committee,  for  their  aid  ,  encouragement,  and 
assistance  in  unraveling  the  sometimes  confusing  jargon  and  mechanisms  involved  in  the 
natural  attenuation  process.  Thanks  are  also  due  to  Lt  Col  Ross  MiUer  of  AFCEE  for  his 
suggesting  the  research  topic  and  for  his  input  to  the  early  development  of  the  research 
project. 

Finally,  I  wi^  to  thank  my  wife,  Kathy,  and  my  children,  Micole,  John,  Aubrey, 
Ruthanne,  and  Shannon  for  the  sacrifice  they  made  by  giving  up  much  of  the  time  they  would 
have  liked  to  spend  with  daddy.  I  am  truly  gratetiil  for  their  patience,  encouragement,  and 
understanding  throughout  this  course  of  study. 


11 


Page 


Preface 

u 

List  of  Figures 

V 

List  of  Tables 

vi 

Abstract 

vii 

I.  Introduction 

1 

General  Issue 

1 

Background 

1 

Specific  Problem 

2 

Objective 

3 

Scope 

4 

Overview 

4 

n.  Literature  Review 

6 

Background 

7 

Definition  of  Terms 

10 

Factors  to  Consider  for  Selection 

24 

BIOPLUME  n™ 

25 

Summary 

28 

m.  Methodology 

30 

Introduction 

30 

Research  Deagn 

30 

Part  I  -  Sensitivity  Analysis 

31 

Part  II  -  Further  Analysis 

37 

IV.  Data  Presentation  and  Discussion 

38 

Introduction 

38 

Part  I 

38 

Part  II 

53 

Some  General  Observations 

59 

iii 


V.  Conclusion  and  Recommemlations 


60 


Introduction  60 

Natural  Attenuation  of  JP-4  in  the  Saturated  Zone  60 

Recommendations  for  Further  Research  61 

Bibliography  63 

Appendix  A  66 

Glossary  of  Terms  67 

Appendix  B  76 

Sample  Input  Data  77 

Appendix  C  78 

Sample  Output  Data  79 

Appendix  D  100 

Data  Input  Formats  101 

Appendix  E  1 10 

Sensitivity  Analysis  111 


iv 


Ua  of.  Figures 

Figure  Page 

1.  Theoretical  and  actual  distributions  of  tetrabydrofiiran  8 

2.  Theoretical  and  actual  distributions  of  phenol  9 

3.  AFCEE  remediation  technology  matrix  12 

4.  Processes  of  dispersion  on  a  microscopic  scale  16 

5.  Relation  between  texture  and  porosity  23 

6.  Simulated  site  and  observation  wells  33 

7.  Changes  in  maximum  concentration  and  distance 

with  changes  in  porosity  40 

8.  Changes  in  maximum  concentration  with  changes 

in  hydraulic  conductivity  43 

9.  Changes  in  distance  with  changes 

in  hydraulic  conductivity  44 

10.  Changes  in  maximum  concentration  and  distance 

with  changes  in  longimdinal  dispersivity  47 

11.  Changes  in  maximum  concentration  and  distance 

with  changes  in  reaeration  coefficient  49 

12.  Changes  in  maximum  concentration  and  distance 

with  changes  in  coefficient  of  anaerobic  biodegradation  52 


V 


List  of  Tables 


Table  Page 

1.  Relative  occurrence  of  benzene,  trichloroethylei»,  and 

tetrachloroethylene  in  water  supply  wells  in  California  10 

2.  Typical  values  of  Porosity  23 

3.  Typical  Hydraulic  Conductivity  Values  24 

4.  Input  data  for  initial  computer  simulation  34 

5.  Data  for  varying  porosity  38 

6.  Data  for  varying  hydraulic  conductivity  41 

7.  Data  for  varying  longitudinal  dispersivity  44 

8.  Data  for  varying  reaeration  coefficient  47 

9.  Data  for  varying  coefficient  of  anaerobic  biodegradation  SO 

10.  Computer  simulations  for  hydraulic  conductivity  of  10*^  ft/s  54 

11.  Computer  simulations  for  hydraulic  conductivity  of  10'^  ft/s  SS 

12.  Computer  simulations  for  hydraulic  conductivity  of  10*^  ft/s  56 

13.  Computer  simulations  for  hydraulic  conductivity  of  10'^  ft/s  57 


VI 


AFIT/GEE/ENV/93S- 12 

Abstract 

This  study  investigated  the  natural  attenuation  mechanisms  and  some  of  the  parameters 
affecting  those  mechanisms  in  the  saturated  zone.  A  literature  search  revealed  numerous 
studies  of  various  attenuation  mechanisms  and  the  associated  parameters.  Much  of  the 
literature  emphasized  biodegradation  as  the  most  promiang  attenuation  mechanism. 
BIOPLUME  II™,  a  fate  and  transport  model,  developed  by  researchers  at  Rice  University  and 
North  Carolina  State  University,  was  used  to  simulate  the  fate  and  transport  of  contaminant 
plume.  The  effects  of  the  model  parameters  were  investigated  by  observing  the  distance  a 
contaminant  plume  was  expected  to  migrate  over  a  fifty  year  period.  The  investigation  was 
limited  by  the  model  which  excludes  chemical  reactions  and  some  physical  and 
physiochemical  reactions.  The  model  simulations  indicated  that  parameters  which  exhibited 
significant  influence  on  natural  attenuation  include  hydraulic  conductivity,  reaeration,  and 
coefficient  of  anaerobic  biodegradation.  These  three  parameters  were  investigated  further  to 
observe  how  they  affected  natural  attenuation  in  concert.  Results  from  the  model  simulations 
indicated  that  the  lower  the  value  of  hydraulic  conductivity  and  the  higher  the  value  of 
aerobic  and  anaerobic  coefficients,  the  more  favorable  are  conditions  for  natural  attenuation. 

A  regression  analysis  was  performed  on  the  data  to  generate  a  mathematical  relationship 
between  site  parameters  and  natural  attenuation. 
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AN  ANALYSIS  OF  SITE  PARAMETERS  AFFECTING 


NATURAL  ATTENUATION  OF  JP-4  IN  THE  SATURATED  ZONE 


1.  InttQdttCtiOP 


General  Issue 

The  mission  of  the  Air  Force  is  "to  defend  the  United  States  through  control  and 
exploitation  of  air  and  space"  (General  Merrill  A.  McPeak,  Air  Force  Chief  of  Staff).  To 
meet  this  mission,  the  Air  Force  must  be  able  to  fly  its  weapon  ^sterns.  One  essential 
requirement  to  getting  these  airplanes  off  the  ground  and  fulfilling  the  mission  is  fuel  for  the 
aircraft.  Therefore,  the  Air  Force  purchases,  transports,  stores,  and  otherwise  handles  vast 
quantities  of  jet  fuel  to  support  their  fundamental  mission.  One  problem  with  possessing  and 
handling  jet  fuel  is  that,  occasionally,  the  fuel  finds  its  way  into  the  soil  by  one  means  or 
another.  This  contamination  can  occur  from  leaking  fuel  storage  tanks  (underground  or  above 
ground),  accidental  fuel  spills,  intentional  fuel  jettisons,  or  some  other  means.  Once  this  fuel 
has  found  its  way  into  the  soil,  it  presents  a  potential  health  hazard  to  the  environment  or  to 
humans.  To  the  extent  this  spent  fuel  poses  a  hazard,  it  must  be  cleaned  up  and  the  affected 
area  restored  to  a  nonhazardous  condition. 

Background 

Several  techniques  have  been  used  in  the  past  to  clean  up  or  remediate  the  sites  where 
jet  fuel  has  made  its  way  into  the  soil.  General  remediation  methods  include  excavation  and 
landfill  or  incineration  and  ex-situ  and  in-situ  chemical  and  biological  treatments  (6:248-254). 
These  remediation  methods  have  met  with  varying  degrees  of  success  and  each  has 
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advantages  and  disadvantages  (31:228-229).  Some  newer  technologies  include  in-situ  heating 
of  soil  by  electrical  energy,  migration  barrier  technology,  air  sparging  of  organic  compounds 
in  soils  and  groundwater,  biological  filters,  and  bioventing  (13:Sessions  1-3).  These  newer 
technologies  exhibit  varying  degrees  of  promise.  Both  the  tried  and  tested  methods  and  the 
newer  technologies  can  be  costly  and  regulatory  compliance  is  often  difficult  to  achieve 
(20:preface).  Because  of  the  magnitude  of  the  problem  and  current  budget  considerations,  it 
is  prudent  to  find  the  least  costly  method  of  remediation.  If  it  could  be  shown  that  no 
remedial  action  was  required  to  eliminate  the  problem,  that  alternative  would  be  a  very 
attractive  option  because  it  would  cost  much  less  to  comply  with  site  closure  criteria.  Money 
saved  could  be  used  to  clean  up  other  sites. 

According  to  Major  Ross  Miller  of  the  Air  Force  Center  for  Environmental 
Excellence  (AFCEE),  natural  attenuation  could  be  a  cost  effective  mechanism  for  remediation 
of  JP-4  at  certain  sites.  He  states  that  from  $1  million  to  $S  million  per  site  could  be  saved 
for  a  total  savings  to  the  Air  Force  of  $500  million  to  $2.5  billion  (23:4). 

Specific  Problem 

The  hydrocarbons  which  reach  the  soil  are  expected  to  attenuate  over  time.  Many 
organic  compounds,  including  hydrocarbons  found  in  petroleum  products,  such  as  JP-4,  are 
broken  down  by  naturally  occurring  bacteria  in  the  soils  (8:482).  The  hazardous  constituents 
of  JP-4  are  the  BTEX  compounds  (13:Session  3),  which  are  also  broken  down  aerobically. 
Under  aerobic  biodegradation,  the  nonhazardous  end  products  of  organic  breakdown  are  water 
and  carbon  dioxide  (1:144).  Hydrocarbons,  contaminating  the  unsaturated  zone,  well  above 
the  water  table,  will  be  subject  to  long-term  aerobic  conditions.  Oxygen  is  all  that  is  needed 
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in  the  unsaturated  zone  to  stimulate  biodegradation  of  the  dissolved  fuel  components.  Due  to 
the  oxygen  content  of  the  air  in  the  unsaturated  zone,  the  hydrocarbons  are  expected  to 
biodegrade  over  time  (8:482).  Biological  clean  up  of  hydrocarbon  contaminated  soils  can  be 
enhanced  by  providing  oxygen  to  the  local  microorganisms  through  forced  aeration  of  the 
soil.  However,  it  is  unnecessary  to  add  oxygen,  other  organisms,  or  nutrients,  in  order  for  the 
attenuation  process  to  proceed  in  the  unsaturated  zone  (23:4). 

Due  to  limited  oxygen  concentrations,  attenuation  in  the  saturated  zone  is  less  certain. 
Because  groundwater  is  the  most  likely  exposure  pathway  for  spilled  or  leaking  JP-4,  it  would 
be  beneficial  to  be  able  to  predict  the  fate  and  transpon  of  the  contaminant  plume  in  the 
saturated  zone.  Knowing  the  expected  behavior  of  the  contaminant  plume  would  allow  one  to 
choose  appropriate  remediation  actions. 

Objective 

Several  models  have  been  developed  to  predict  degradation  of  various  contaminants  in 
certain  environments.  There  are  ongoing  studies  to  validate,  calibrate,  and  improve  these 
models.  The  model  chosen  for  use  in  this  study  is  the  BIOPLUME  II  model. 

The  focus  of  this  research  is  to  investigate,  through  the  use  of  the  BIOPLUME  II 
model,  the  conditions  and  parameters  which  favor  natural  attenuation,  to  examine  the  plume 
characteristics  for  these  conditions  and  parameters,  and  to  theoretically  estimate  whether  the 
contaminant  concentrations  will  fall  within  established  risk  and  contamination  limits  before 
reaching  a  potential  exposure  point  and  becoming  a  hazard  to  human  health  or  to  the 
environment.  If  risk  and  contaminant  levels  fall  within  established  guidelines,  the  natural 
attenuation  alternative  may  be  the  appropriate  remediation  elective. 
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This  research  project  proposes  to  use  BIOPLUME  II,  software  developed  by  Rice 
University,  to  assess  the  site  conditions  and  parameters  which  enhance  or  optimize  natural 
attenuation  of  JP-4.  Theoretical  data  will  be  input  to  the  model  to  analyze  how  the  value  of 
various  parameters  affect  the  contaminant  plume.  A  sensitivity  analysis  will  be  accomplished 
to  explore  critical  ranges  and  variations  in  parameters  and  to  explore  which  parameters  may 
be  the  most  critical.  General  site  conditions  will  be  proposed  for  which  natural  attenuation 
should  be  considered  the  primary  remediation  option. 

Scope 

Natural  attenuation  includes  several  mechanisms  by  which  contaminant  concentrations 
are  reduced.  These  mechanisms  include  volatilization,  sorption,  dilution,  dispersion, 
advection,  hydrolysis,  biodegradation,  and  others.  These  processes  occur  in  distinct  zones, 
three  of  which  are  the  atmosphere,  the  vadose  or  unsaturated  zone,  and  the  saturated  zone. 
This  research  will  focus  on  natural  attenuation  in  the  saturated  zone  which  includes  natural 
biotransformation  of  hazardous  organic  compounds  in  the  subsurface  (37:1).  The 
investigation  will  be  limited  to  the  analysis  of  results  generated  by  the  BIOPLUME  II™ 
software. 

Overview 

This  chapter  has  given  a  brief  background  to  set  the  stage  for  the  subsequent 
discussion  of  natural  attenuation  and  associated  parameters.  These  parameters  will  be 
examined  with  the  aid  of  BIOPLUME  II™,  a  two-dimensional  model  which  will  simulate  the 
natural  attenuation  process  in  the  saturated  zone.  Chapter  2  of  the  thesis  consists  of  a 
discussion  of  the  relevant  literature  on  this  topic.  Chapter  3  includes  the  research  design  and 
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the  methodology  employed  for  this  research  effon.  Chiq|)ter  4  contains  the  data  description 
and  analysis  and  chapter  S  presents  discussion  of  findings,  conclusion,  and  recommendations. 
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In  the  last  decade  or  two,  the  Resource  Conservation  and  Recovery  Act  (RCRA)  and 
the  Comprehensive  Environmental  Response,  Compensation,  and  Liability  Act  (CERCLA) 
have  prompted  a  £ren2y  of  activity  in  the  environmental  restoration  arena.  Because  of  the 
increased  attention  given  to  the  environment,  a  great  deal  of  research  has  been  done  in  the 
Eeld  of  hazardous  waste  remediation.  Due  to  the  increasing  costs  of  environmental 
restoration,  researchers  have  been  looking  for  effective  remediation  technologies  which  are 
also  cost  effective.  The  search  has  led  to  few  promising  technologies  which  are  effective  yet 
inexpensive.  One  of  the  least  costly  methods  of  restoration  is  natural  attenuation.  Even 
where  naturally  occurring  mechanisms  are  at  work  degrading  hazardous  substances  into 
harmless  byproducts,  one  must  verify  that  this  activity  is  taking  place  and  that  the  natural 
process  will  be  sufficient  to  bring  the  hazard  within  acceptable  limits.  This  verification  is 
accomplished  by  monitoring  the  contaminant  plume  at  a  hazardous  waste  site  through  the  use 
of  monitoring  wells.  Even  with  extensive  monitoring  requirements,  cost  savings  can  be 
significant.  Natural  attenuation  has  become  a  very  promising  remediation  option  due  to  its 
low  cost.  Wherever  possible,  the  Air  Force  would  like  to  take  advantage  of  this  natural 
degradation  process  because  of  the  attractive  low  cost  advantage  (23  ;4). 

This  chapter  will  take  a  closer  look  at  the  natural  attenuation  alternative.  The  chapter 
will  begin  with  some  background  followed  by  a  deHnition  of  some  of  the  terms  associated 
with  natural  attenuation  and  a  brief  description  of  each  of  the  principle  mechanisms 
comprising  the  natural  attenuation  process.  Following  the  discussion  of  natural  attenuation 


mechanisms,  the  soil  parameters  or  aquifer  parameters  affecting  these  natural  mechanisms  will 
be  addressed.  The  discussion  will  be  limited  to  activity  in  the  saturated  zone.  Next,  factors 
to  be  considered  in  selecting  natural  attenuation  as  a  remediation  option  will  be  reviewed.  A 
brief  overview  of  BIOPLUME  II™,  the  model  chosen  to  analyze  the  natural  attenuation 
parameters,  will  be  presented  followed  by  a  final  summary. 

BacisgQUBd 

Early  studies  of  contaminant  transport  indicated  that  contaminant  plumes  did  not 
migrate  as  fast  and  as  far  as  initially  predicted  (27:208).  Figures  1  and  2  show  the 
differences  between  predicted  [theoretical]  contaminant  plumes  and  actual  [measured] 
contaminant  plumes  for  two  speciHc  compounds  at  one  site.  Observations  like  this  led  to  the 
discovery  of  various  naturally  occurring  attenuation  mechanisms.  More  recent  studies  have 
been  conducted  to  better  understand  the  natural  attenuation  mechanisms.  An  aerobic 
biodegradation  study  of  water  supply  wells  in  California,  conducted  by  Hadley  and  Armstrong 
in  1991,  compared  the  presence  of  benzene,  an  eaaly  biodegraded  contaminant,  with  two 
solvents,  trichloroethylene  (TCE)  and  tetrachloroethylene  (PCE),  which  are  not  as  easily 
degraded.  The  study  indicated  that  far  less  benzene  enters  the  subsurface  groundwater  than 
expected.  Table  1  shows  the  relative  occurrence  of  benzene  as  compared  to  the  other  two 
solvents  in  the  California  wells.  The  results  from  these  studies  suggest  that  naturally 
occurring  mechanisms  are  at  work  removing  contaminants  from  groundwater  under  certain 
conditions  (38:2). 

Much  has  been  written  concerning  the  possible  mechanisms,  including  volatilization, 
sorption,  chemical  and  biological  degradation,  and  others,  that  come  into  play  in  the 
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1.  Theoretical  and  actual  di^butions  of  phenol  (27:213) 


Table  1.  Relative  occurrence  of  benzene,  trichloroethylene,  and  tetrachloroethylene  in  water 
supply  wells  in  California.  (38:2) 


Compound 

Number  of 
Wells 

Median 

Concentration 

Maximum 

Concentration 

Minimum 
Reported  Cone. 

(mg/liter) 

(mg/liter) 

(mg/liter) 

Benzene 

9 

0.2 

1.1 

0.1 

TCE 

188 

3.2 

538 

0.1 

PCE 

199 

1.9 

166 

0.1 

attenuation  of  chemical  contaminants.  Research  has  shown  that  under  varying  site  conditions, 
great  differences  occur  between  attenuation  efficiencies  observed  in  the  laboratory  and  those 
encountered  in  the  Held  (24:830).  This  variation  may  be  attributable  to  the  complexity  of  the 
field  environment  and  the  many  parameters  that  affect  expected  attenuation  efficiencies 
(20:205). 

Definition  of  Terms 

In  the  following  discussion,  the  principle  mechanisms  associated  with  natural 
attenuation  will  be  defined  and  discussed.  Many  of  the  terms  associated  with  natural 
attenuation  not  defined  here  may  be  found  in  a  glossary  of  often-encountered  terms  in 
Appendix  A. 

Attenuation 

Attenuation  of  chemical  constituents  in  soil  and  rock  is  a  well  recognized 
phenomenon.  The  concept  has  been  and  is  currently  being  used  by  state  and  federal 
regulators  for  the  following  actions  (20:215): 
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1)  establishing  general  standards 

2)  determining  cleanup  criteria  for  polluted  sites  on  a  case  by  case  basis 

3)  providing  variances  for  the  disposal  of  waste  into  landfills 

4)  regulating  the  disposal  of  treated  waste  water  sludge  to  land 

5)  regulating  ground  water  recharge  projects. 

"Attenuation  is  defined  as  the  reduction  in  quantity  or  toxicity  of  a  chemical  constituent  in 
liquids  that  may  be  released  to  the  environment  (20:205)."  Attenuation  can  be  accelerated 
through  various  remediation  technologies;  however,  this  study  is  concerned  with  the  natural, 
rather  than  anthropogenic,  processes. 

Natural  Attenuation 

Namral  attenuation  is  the  naturally  occurring  attenuation  of  a  chemical  constituent 
without  artificial  stimulation  or  augmentation.  Natural  attenuation  occurs  because  of  several 
mechanisms.  These  mechanisms  have  been  delineated  into  physical,  physiochemical, 
chemical,  and  biological  categories.  These  mechanisms  operate  along  the  contaminant  plume 
migration  pathway  and  begin  almost  immediately  when  a  contaminant  is  released  to  the 
environment  (20:217). 

Natural  attenuation  is  the  preferred  method  of  remediation  in  the  Air  Force  as  shown 
in  Figure  3.  This  matrix  presents  a  hierarchy  of  preferred  remediation  alternatives.  As 
shown,  natural  attenuation  is  the  preferred  method  in  all  situations  except  vapor  treatment  for 
which  natural  remediation  is  not  applicable. 


.  AFCEE  remediation  technology  matrix  (24:4) 


Physical  Mechanisms 

The  phyacal  mechanisms  to  be  discussed  include  advection,  filtration,  volatilization, 
dispersion,  and  diffusion.  These  physical  mechanisms  attenuate  a  contaminant  primarily  by 
reducing  the  concentration  of  contaminant  rather  than  reducing  the  toxicity  of  the 
contaminant.  Frequently,  the  physical  mechanisms  simply  transfer  the  contaminant  from  one 
medium  to  another  (21:5). 

Advection 

In  permeable  aquifers,  such  as,  sand  and  gravel  aquifers,  advection  is  the  primary 
mechanism  affecting  contaminant  migration  (9:159).  Advection  is  a  result  of  the  bulk  motion 
of  flowing  groundwater  (29:1-6).  The  flowing  groundwater  transports  a  contaminant  through 
the  substrate  from  a  region  of  high  hydraulic  he^  to  a  region  of  lower  hydraulic  head 
(9:159).  The  difference  in  hydraulic  head  causes  a  pressure  differential  termed  the  hydraulic 
gradient.  The  velocity  at  which  groundwater  flows  through  a  substrate  is  a  function  of  the 
hydraulic  gradient  and  the  substrate's  inherent  capability  to  transport  groundwater.  This 
inherent  capability  is  a  function  of  the  soil  or  aquifer  characteristics  which  will  be  discussed 
later.  Calabresi  and  Kostecki  report  that  "Groundwater  velocities  in  uniform  sand  and  gravel 
aquifers,  under  natural  gradient  conditions,  are  typically  between  10  and  100  m/yr,  with  a 
potential  range  between  1  and  1000  m/yr  (9:159)." 

Filtration 

Filtration  is  defined  as  the  simple  obstruction  of  particulates  migrating  through  a 
permeable  soil.  The  particulates  may  be  insoluble  minerals  or  chemical  precipitates.  Physical 
filtration  is  the  process  of  suspended  solids  being  captured  in  the  soil  matrix  because  they  are 
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too  large  to  flow  through  the  pore  spaces  (20:222)  Filtration  attenuates  contaminants  by 
obstructing  their  movement  through  the  substrate.  A  change  in  groundwater  pH  or  other 
environmental  conditions  may  cause  a  chemical  precipitate  or  insoluble  mineral  to  dissolve 
and  be  released  to  the  groundwater.  For  this  reason,  filtration  should  not  be  relied  on  as  an 
attenuation  mechanism. 

Volatilization 

Volatilization  is  defined  as  the  process  whereby  a  contaminant  is  transferred  from  the 
soil  or  water  phase  to  air  or  vapor  phase.  Volatilization  is  a  form  of  diffusion,  that  is, 
molecular  movement  of  chemicals  or  ions  from  a  high  concentration  region  to  a  low 
concentration  region.  Volatilization  is  an  extremely  important  attenuation  mechanism  for 
many  organic  compounds;  however,  rates  of  volatilization  can  vary  over  a  large  range  (9:160). 

In  the  saturated  zone,  there  is  no  air  in  the  soil  pores,  therefore,  only  the  surface  of  the 
groundwater  table  or  capillary  fringe  is  in  contact  with  air.  Also,  the  aquifer  is  not  a  well- 
mixed  system.  Lyman  reports  that  these  two  conditions  eliminate  volatilization  from  being  a 
significant  attenuation  mechanism  in  the  saturated  zone  (23:209). 

Dispersion 

Dispersion  refers  to  the  spreading  out  of  a  dissolved  contaminant  as  it  is  transported 
by  the  groundwater  through  the  soil  matrix.  Dispersion  occurs  because  of  two  fundamental 
processes,  molecular  diffusion  in  solution  due  to  a  concentration  gradient  and  mechanical 
mixing  (or  mechanical  dispersion)  caused  by  the  groundwater  flow  through  the  soil.  The 
molecular  diffusion  process  is  a  result  of  the  contaminant  seeking  chemical  equilibrium.  The 
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overall  result  is  a  transfer  of  solutes  from  a  high  concentration  zone  to  a  lower  concentration 
zone  (9:  159). 

Molecular  diffusion  is  a  process,  whereby  ionic  or  molecular  constituents  move  under 
the  influence  of  their  kinetic  activity  causing  a  net  directional  movement  opposite  the 
direction  of  their  concentration  gradient.  The  mass  of  a  contaminant  is  diffused  as  it  passes 
through  a  specified  cross  section  of  a  given  medium  per  unit  of  time.  According  to  Pick's 
first  law,  the  diffusion  of  a  contaminant  is  proportional  to  the  concentration  gradient. 

The  rate  at  which  a  chemical  contaminant  is  diffused  is  a  function  of  the  chemical 
properties  of  both  the  pollutant  and  the  soil  through  which  it  passes.  In  the  saturated  zone, 
the  effect  of  mechanical  dispersion  is  generally  much  greater  than  that  of  molecular  diffusion, 
because  mechanical  dispersion  is  proportional  to  the  groundwater  flow  velocity.  Therefore, 
diffusion  is  often  neglected  in  the  saturated  zone.  Diffusion  is  generally  not  significant  unless 
the  aquifer  has  a  very  low  flow  rate.  Tucker  and  Nelken  suggest  that  molecular  diffusion  can 
be  ignored  at  groundwater  velocities  exceeding  .002  cm/sec  (9:159-160). 

Dispersion  can  have  a  significant  effect  on  natural  attenuation  in  the  saturated  zone 
because  it  exposes  the  chemical  contaminant  to  a  larger  volume  of  soil  and  pore  water  at  a 
reduced  concentration,  increasing  the  potential  for  attenuation  by  sorption  or  degradation  by 
chemical  reactions  and  microbial  activity  (20:223).  Figure  4  illustrates  the  dispersion 
processes  on  a  molecular  scale. 
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Figure  4;  Processes  of  Dispersion  on  a  Microscopic  Scale  (30: 1  -  8) 


Physiochemical  Mechanisms 

The  mechanisms  included  in  this  category  are  neither  purely  physical  nor  purely 
chemical  but  depend  on  both.  The  primary  mechanism  included  in  this  category  is  sorption. 
Sorption 

Sorpdon  consists  of  absorption,  adsorption,  and  desorption.  Absorption,  or  solid  state 
diffusion,  is  the  accumulation  of  a  contaminant  within  the  soil  particle  pores.  Absorption 
results  in  a  decrease  in  the  amount  of  an  element  available  for  transport  in  the  groundwater, 
effectively  reducing  the  concentration  of  contaminant  in  the  groundwater.  Adsorption  is  the 
accumulation  of  the  pollutant  at  the  surface  of  the  soil  particle.  The  bonding  which  causes 
this  accumulation  can  be  either  physical  or  chemical.  The  chemical  bonding  is  also  referred 
to  as  chemisorption  and  is  usually  irreversible,  rendering  the  chemical  immobile.  Adsorption 
can  decrease  the  concentration  of  an  element  in  the  groundwater.  Desorption  is  transfer  of 
the  contaminant  from  the  soil  particle  pores  to  the  soil  particle  surface  or  from  the  soil 
particle  surface  back  into  the  groundwater  (9:160-163,  20:228). 
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One  significant  effect  sorption  has  in  the  attenuation  of  contaminants  is  that  it  can 
drastically  retard  the  migration  of  pollutants  in  soils.  Sorption  and  subsequent  migration 
retardation  are  dependent  upon  the  type  of  chemical  and  the  type  of  soil  through  which  a 
contaminant  is  moving  (9:163). 

Sorption,  like  filtration,  should  not  be  relied  upon  as  an  attenuation  mechanism. 
Because  of  the  retardation  effect,  sorption  can  aid  in  controlling  the  migration  of  a 
contaminant  plume;  however,  post  remediation  monitoring  for  certain  remediation  techniques 
has  shown  that  sorption  may  be  reversible  after  a  site  has  been  restored  (28:42-44).  The  low 
contaminant  concentrations,  resulting  from  remediation,  cause  subsequent  desorption  of  the 
contaminant  from  the  soil  particles  back  into  the  groundwater  resulting  in  elevated 
contaminant  concentrations. 

Chemical  Mechanisms 

Chemical  mechanisms  can  result  in  transformation  of  the  contaminant  from  a 
hazardous  chemical  species  to  a  nonhazardous  species  (20:233).  These  mechanisms  are 
considered  equilibrium  processes  and  may  be  reversed  if  the  physiochemical  conditions  of  the 
pore  water  change  (20:222).  The  reversibility  of  the  chemical  reactions  is  one  reason  certain 
chemical  mechanisms  should  not  be  relied  upon  as  attenuation  mechanisms.  Chemical 
mechanisms  include  precipitation,  hydrolysis,  complexadon,  oxidation,  and  reduction. 
Precipitation 

Precipitation  occurs  when  the  contaminant  exceeds  its  solubility  limits  and  forms  an 
insoluble  solid  compound  which  becomes  suspended  in  solution.  The  solid  particles  are  then 
subject  to  the  filtering  mechanism  discussed  earlier.  Filtration  actually  "captures"  the 
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constituent  in  the  soil  matrix,  preventing  it  hrom  migrating.  These  suspended  particles  can 
remain  insoluble  indefinitely;  however,  a  change  in  environmental  conditions,  such  as  redox 
potential,  pH,  or  ion  concentration,  may  cause  the  precipitate  to  dissolve  and  be  released  back 
into  the  groundwater.  Because  of  the  reversibility  of  these  chemical  reactions,  precipitation 
should  not  be  relied  upon  as  an  attenuation  mechanism  (20:231-232). 

Hydrolysis 

A  hydrolysis  reaction  results  in  the  transformation  of  contaminants  from  one  chemical 
species  to  another.  Hydrolyas  reactions  primarily  affect  organic  compounds.  Organic 
compounds  are  transformed  as  they  react  with  water,  forming  new  compounds.  The  process 
typically  consists  of  the  breaking  of  the  carbon-halogen  bond  in  the  original  molecule  and  the 
formation  of  a  new  carbon-oxygen  bond  (9:166-168). 

Depending  on  the  chemical  compound  and  the  environmental  conditions,  the  rate  of 
hydrolysis  can  vary  widely.  Hydrolysis  is  sit'?'  dependent  and  may  or  may  not  be  a  significant 
attenuation  mechanism  in  the  saturated  zone.  Fourteen  orders  of  magnitude  and  associated 
half-lives  as  low  as  a  few  seconds  to  as  high  as  10‘  years  can  be  expected  (9:167).  Because 
of  this  variability,  hydrolysis  is  often  excluded  from  models  or,  if  included,  it  is  often  not 
corrected  for  environmental  conditions  (9:167).  In  either  case,  model  results  should  be 
evaluated  cautiously.  The  BIOPLUME  II™  model  used  in  this  study  excludes  hydrolysis 
from  consideration,  therefore,  for  some  conditions,  the  model  may  yield  conservative  results. 
Complexation 

Complexation,  or  chelation,  occurs  when  a  ch«;mical  complex  is  formed  by  the  union 
of  a  metal  ion  with  nonmetallic  ions.  The  resultant  metal-ligand  complex  will  generally 
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prevent  the  metal  from  further  chemical  reactions  (9:169).  The  complex  may  also  be  less 
toxic  than  either  of  the  complexed  chemicals  (^:234).  Because  it  is  primarily  an  attenuation 
mechanism  for  metals,  complexation  is  probably  not  a  significant  attenuation  mechanism  for 
jet  fuel  in  the  saturated  zone.  Due  to  limited  understanding  and  available  information  about 
the  process,  complexation  reactions  have  not  been  modeled  (9:169). 

Oxidation-Reduction 

For  some  organic  compounds,  chemical  oxidation  is  an  important  degradation  process. 
The  process  is  dependent  upon  environmental  conditions  and  usually  involves  reactions 
between  the  organic  compound  and  free  radicals  already  in  solution  (9:166).  Oxidation 
occurs  when  an  element  undergoes  a  change  in  valence  state,  from  a  lower  valence  state  to  a 
higher  valence  state,  by  the  loss  of  valence  electrons.  The  change  in  valence  often  transforms 
the  element  from  a  solid  state  to  a  soluble  ion.  Reduction  occurs  when  the  chemical 
contaminant  gains  valence  electrons.  Either  process  may  result  in  reduced  toxicity  of  a 
chemical  contaminant.  Reduction  in  toxicity  occurs  when  the  valence  electron  state  of  a 
contaminant  causes  the  chemical  species  to  become  less  toxic  than  the  original  chemical 
species.  This  electron  valence  change  can  also  alter  the  mobility  of  the  contaminant  (20:234). 
Biological  Mechanisms 

By  far,  the  majority  of  literature  pertaining  to  natural  attenuation  emphasizes 
bioattenuation.  Several  studies  have  shown  that  naturally  occurring  microorganisms  can 
degrade  or  oxidize  hazardous  substances  such  as  jet  fuel  into  harmless  substances  such  as 
water  and  carbon  dioxide  (13:Sessions  2-  4,  37:380,  9:115).  This  process,  known  as  natural 
biodegradation,  is  fast  becoming  a  preferred  remediation  alternative.  One  author  states  that 
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"Bioremediation,  or  enhanced  microbiological  treatinent,  of  environments  contaminated  with  a 
variety  of  organic  and  inorganic  compounds  is  one  of  the  most  effective  innovative 
technologies  to  come  along  this  century"  (21:preface). 

In  reality,  bioremediation  is  merely  a  new  application  of  a  technology  used  in  ancient 
times  for  wastewater  treatment.  "Romans  and  others  built  intricate  networks  of  sewers  as 
early  as  600  B.C.  for  collecting  wastewater  which  underwent  subsequent  biological  treatment" 
(21:1).  Natural  biodegradation,  however,  existed  long  before  these  manmade  uses  of 
microorganisms. 

Biodegradation 

Biodegradation  in  the  subsurface  includes  three  patterns  of  biotransformation:  aerobic 
respiration,  anaerobic  respiration,  and  fermentation  (37:368,  38:1,9:115).  The  best  understood 
and  most  developed  bioremediation  technologies  rely  on  aerobic  microbial  activity.  Aerobic 
degradation  is  considered  the  most  attractive  degradation  pattern  due  to  its  rapid  rate  of 
transformation  and  the  absence  of  noxious  by-products  which  may  result  from  anaerobic 
microbial  degradation  (37:368).  In  aerobic  respiration,  microorganisms  consume  oxygen 
while  decomposing  organic  hydrocarbons  (26:73).  The  basic  requirements  for  aerobic 
biodegradation  include:  the  abundance  of  local  microorganisms;  the  availability  of  an  energy 
source,  a  carbon  source,  oxygen  and  nutrients;  and  the  presence  of  acceptable  environmental 
conditions  (6:Rifai).  Dissolved  oxygen  concentration  in  the  groundwater  is  usually  the 
limiting  factor  for  aerobic  respiration  (23:296).  Natural  aerobic  biodegradation  can  account  for 
significant  mass  loss  of  dissolved  contaminants  in  a  plume  (13:Rifai). 
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Anaerobic  decomposition  of  hazardous  chemicals  occurs  in  the  absence  of  oxygen.  In 
an  oxygen  depleted  environment,  other  electron  acceptors  include  nitrate,  sulfate,  carbonate, 
and  iron  in  (9:115).  The  effectiveness  of  anaerobic  degradation  is  somewhat  more 
controversial  than  aerobic  degradation.  Calabresi  and  Kostecki  state  that  "convincing  proof  of 
significant  anaerobic  hydrocarbon  biodegradation  is  still  outstanding.  Sulfates  are  a  potential 
electron  acceptor,  but  are  not  abundant  in  soils.  Nitrate  is  not  energetically  favorable  for  this 
purpose  in  soils  (8:224)."  Although  this  process  is  not  yet  well  understood,  recent  studies 
show  that  certain  chemical  compounds  are  biodegradable  in  groundwater  under  anaerobic 
conditions.  Anaerobic  degradation  is  a  slower  process  than  aerobic  degradation,  sometimes 
requiring  months  to  years  for  the  process  to  begin  (38:4). 

In  the  fermentation  process,  the  organic  compound  serves  as  an  electron  acceptor, 
instead  of  the  more  common  roles  of  electron  donor  or  carbon  source.  In  essence,  the 
chemical  contaminant  becomes  something  for  the  microorganism  to  breath  rather  than  a  food 
source.  In  this  process,  the  degradation  rate  is  limited  by  the  electron  acceptor  demand  from 
other  organic  compounds  and  by  competition  from  other  more  conventional  electron  acceptors 
(38:1-4). 

Aquifer  Parameters 

The  substrate  properties  and  characteristics  which  are  considered  important  to  physical 
transport  in  the  saturated  zone  include  composition  of  the  soil,  storage  capacity  of  the  soil, 
porosity,  hydraulic  conductivity,  and  dispersivity  parameters. 
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Porosity 

Porosity  is  defined  as  the  ratio  of  the  volume  of  the  void  space  between  soil  particles 
to  the  total  volume  by  the  relationship: 


n  = 
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voids 


V 
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Figure  S  illustrates  the  relationship  between  the  soil  texture  and  porosity. 

Typical  values  for  porosity  are  given  in  Table  2.  A  common  misconception  is  that  the  more 
porous  the  soil,  the  greater  the  ability  to  transport  groundwater.  However,  one  can  see  that 
porosity  varies  over  a  wide  range  for  different  soil  types,  such  as  well  sorted  sand  and  gravel 
or  clay.  Clay  soils  can  have  a  greater  porosity  than  sand  and  gravel,  yet  are  less  permeable, 
that  is,  less  able  to  transport  groundwater. 

Dispersivity 

Dispersion  has  been  discussed  previously  as  an  attenuation  mechanism.  The  mechanism, 
however,  is  a  function  of  soil  structure  and  composition.  The  spreading  of  contaminant  in  the 
same  direction  as  groundwater  flow  is  termed  longitudinal  dispersion.  Dispersion 
perpendicular  to  the  direction  of  flow  is  termed  transverse  dispersion  (29:1-8). 
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Figure  5.  Relation  between  texture  and  porosity,  (a)  Well-sorted  sedimentary  deposit  having 
high  porosity;  (b)  poorly  sorted  sedimentary  deposit  having  high  porosity;  (c)  well-sorted 
sedimentary  deposit  consisting  of  pebbles  that  are  themselves  porous,  so  that  the  deposit  as  a 
whole  has  a  very  high  porosity;  (d)  well-sorted  sedimentary  deposit  whose  porosity  has  been 
diminished  by  the  deposition  of  mineral  matter  in  the  interstices;  (e)  rock  rendered  porous  by 
solution;  (f)  rock  rendered  porous  by  fracturing.  (12:25) 


Table  2.  Typical  values  of  Porosity  (30:1  -  3) 


Well  sorted  sand  or  gravel 

25  -  50%  1 

Sand  and  gravel,  mixed 

20  -  35% 

Glacial  till 

10  -  20% 

SUt 

35  -  50% 

Clay 

33  -  60% 
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Hydraulic  Conductivity 


Permeability  is  the  term  which  refers  to  the  capacity  of  soil  to  transmit  a  fluid  and  is 
related  to  hydraulic  conductivity.  Hydraulic  conductivity  refers  to  the  rate  at  which  an  aquifer 
conveys  water.  Hydraulic  conductivity  values  are  usually  obtained  from  testing  at  a  site. 
Typical  values  for  various  soil  types  are  included  in  Table  3. 


Table  3.  Typical  Hydraulic  Conductivity  Values  (cm/sec)  (30:1-7) 


Clay 

10  ’  -  10^ 

Silt,sandy  silts,  clayey  sands,  tills 

10“"  -  lO"* 

Silty  sands,  fine  sands.  Well-sorted  sands. 
Glacial  outwash 

10^  -  10' 

Well-sorted  gravel 

10^  -  1 

Factors  to  Consider  for  Selection 

Factors  that  affect  the  natural  attenuation  process  in  the  saturated  zone  need  to  be 
understood.  General  site  parameters  that  must  be  investigated  when  considering  the 
applicability  of  natural  attenuation  include  (8:178): 

-  Soil  microbiology:  First  of  all,  microorganisms  which  have  the  ability  to 
degrade  petroleum  products  must  exist  in  the  zone  of  contamination. 

-  Soil  chemistry:  In  order  to  maintain  a  healthy  population  of 
microorganisms,  not  only  must  oxygen  and  nitrogen  supplies  be 
sufficient  but  other  soil  nutrients  must  also  be  in  adequate  supply. 
Furthermore,  conditions  which  preclude  a  healthy  environment  must  not 
exist,  such  as  toxic  salts  or  heavy  metals  in  the  soil. 
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-  Soil  physics:  The  soil  must  be  porous  and  permeable  enough  for  oxygen 
transfer  to  take  place.  Movement  of  nutrients,  including  oxygen  and 
nitrogen,  into  the  area  of  contamination  must  proceed,  as  well  as 
transfer  of  carbon  dioxide  away  from  the  contaminated  area. 

-  Soil  morphology:  Understanding  of  the  subsurface  geology  is  necessary 
for  an  effective  remediation  design. 

-  Hydrogeology:  As  much  information  as  possible  should  be  gathered 
concerning  aquifer  characteristics:  depth  to  the  water  table,  hydraulic 
gradient,  direction  of  flow,  and  initial  concentrations  of  nutrients  and 
contaminants.  The  smaller  the  contaminant  source,  the  more  likely  it  is 
that  natural  attenuation  will  be  sufficient  to  degrade  the  contaminant 
plume. 

Other  parameters  that  help  assess  whether  or  not  natural  attenuation  ought  to  be 
considered  as  a  remediation  option  include  the  following  (17:117-118): 

-  suitability  of  groundwater  for  consumption 

-  potential  for  exposure,  i.e.  distance  to  exposure  point 

-  feasibility  of  active  restoration  techniques 

-  projected  future  demand  for  groundwater 

-  proximity  to  a  surface  water  discharge  area  with  subsequent  dilution  to 
environmentally  safe  levels 

-  attainability  of  acceptable  contaminant  and  risk  levels. 


Noonan  includes  sociopolitical  considerations  as  another  factor  that  ought  to  be 

considered  in  the  remediation  design  and  argues  that: 

Because  in  situ  treatment  [natural  attenuation]  is  accomplished 
largely  underground,  little,  if  any,  evidence  of  activity  may  be 
discernible  by  the  general  public.  This  lack  of  observable  activity 
may  lead  to  a  public  perception  of  no  action.  The  public  in  general 
may  be  more  prone  to  respond  positively  to  corrective  actions  that 
manifest  significant  levels  of  activity.  (26:83) 


In  order  to  derive  maximum  benefit  from  the  namral  attenuation  process,  it  would  be 
advantageous  to  be  able  to  predict  which  site  conditions  favor  natural  attenuation,  thereby 
indicating  which  sites  would  be  candidates  for  natural  attenuation  as  the  preferred 
remediation  method.  In  this  study,  it  is  proposed  to  examine  the  parameters  involved  in 
natural  attenuation  for  the  purpose  of  predicting  general  site  conditions  which  favor  the 
natural  attenuation  process.  Investigation  of  these  parameters  and  their  effect  on  the  natural 
attenuation  process  will  be  accomplished  using  the  BIOPLUME  II™  model. 

BIQRLUMEieL 

The  purpose  of  this  research  is  threefold:  1)  to  investigate,  through  the  use  of  the 
BIOPLUME  II™  model,  the  conditions  and  parameters  which  favor  natural  attenuation,  2)  to 
assess  the  degradation  rate  for  these  conditions  and  parameters,  and  3)  to  theoretically 
estimate  whether  this  degradation  rate  causes  the  contaminant  concentrations  to  fall  within 
establitiied  risk  and  contamination  limits  before  reaching  a  potential  exposure  point  and 
becoming  a  hazard  to  human  health  or  the  environment.  If  risk  and  contaminant  levels  fall 
within  established  guidelines,  the  no-action  or  natural  attenuation  alternative  may  be  the 
appropriate  remediation  elective. 

BIOPLUME  n™  is  an  extremely  versatile,  two-dimensional  computer  model  which 
simulates  the  migration  and  attenuation  of  dissolved  hydrocarbon  plumes.  The  model 
includes  aerobic  and  anaerobic  biodegradation,  sorption,  and  dispersion  of  the  contaminant 
plume.  The  model  is  also  able  to  simulate  reaeration  and  radioactive  decay.  Reaeration  and 
the  anaerobic  biodegradation  process  are  modeled  as  first  order  reactions  (28:v). 
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An  input  data  set  is  presented  in  Appendix  B.  This  data  set  is  simplified  and  input 
data  can  be  varied  to  consider  parameters  of  particular  interest.  The  output  data  for  the  model 
is  much  more  extensive.  An  example  is  provided  as  Appendix  C.  The  output  file  includes  all 
the  input  data,  followed  by  the  simulation  results.  The  length  of  the  output  is  dependent  on 
the  input  data  and  can  be  very  lengthy.  Specific  input  parameters  are  included  as  Appendix 
D. 

The  model  calculates  changes  in  the  contaminant  and  oxygen  concentrations  for  a 
given  time  period  by  solving  the  solute  transport  equation  for  the  contaminant  and  oxygen 
plumes  and  then  superimposing  the  results.  The  model  assumes  an  instantaneous  reaction 
between  oxygen  and  hydrocarbon  to  simulate  the  biodegradation  processes  (28  :v).  This 
assumption  is  evident  in  viewing  the  output.  Even  though  contaminant  effects  are  evident  in 
a  cell,  by  observing  the  oxygen  reduction  taking  place,  contaminant  concentrations  are  not 
displayed  as  an  increase  or  accumulation  in  a  cell  until  oxygen  concentrations  are  zero. 

BIOPLUME  II™  includes  three  sources  of  oxygen  into  an  aquifer:  initial  dissolved 
oxygen  in  the  uncontaminated  aquifer,  natural  recharge  of  oxygen  across  a  vertical  cross- 
section,  and  vertical  exchange  of  oxygen  from  the  unsaturated  zone  (reaeration)  (28  :v).  These 
three  sources  can  be  used  to  provide  oxygen  to  drive  natural  biodegradation  processes  to 
attenuate  a  contaminant  plume.  In  the  event  one  desires  to  investigate  "active"  bioremediation 
techniques,  injection  wells  can  be  specified  as  oxygen  sources  into  a  contaminated  aquifer  as 
well(28:v). 

The  BIOPLUME  II™  model  has  been  used  effectively  at  several  ates.  The  model 
provided  an  adequate  simulation  of  a  plume  of  contamination  from  a  forty  year  old  creosote 
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waste  lagoon.  The  plume  had  undergone  .  gnificant  attenuation  over  twenty  years,  migrating 
only  400  meters.  BIOPLUME  simulations  indicated  that  the  contaminant  plume  had 
reacned  its  maximum  extent  (38:3). 

In  another  study,  the  model  was  applied  to  an  alkylbenzene  plume  emanating  from  a 
flare  pit  of  a  natural  gas  facility.  Again,  the  model  adequately  predicted  contaminant 
concentrations  at  the  site  (38:3) 

A  sensitivity  analysis  conducted  by  Hanadi  Rifai,  author  of  the  model,  is  provided  as 
Appendix  E.  The  sensitivity  analysis  relates  the  parameters  in  the  model  to  the  natural 
attenuation  process.  It  was  concluded  that  naturally  attenuated  contaminant  mass  in  the  model 
is  most  sensitive  to  hydraulic  conductivity,  the  coefficient  of  reaeration,  and  the  coefficient  of 
anaerobic  decay  (28:2-8).  In  this  study,  to  add  further  insight  to  plume  behavior,  these  three 
parameters  were  investigated  more  fully.  Chapter  three  provides  a  detailed  methodology  of 
this  investigation. 

Summary 

Due  to  various  mechanisms,  natural  attenuation  begins  almost  immediately  upon 
release  of  a  contaminant  to  the  environment.  (6:sessions  1-4).  In  general,  physical 
mechanisms  have  a  lower  attenuation  potential;  however,  they  can  be  important  in  support  of 
physiochemical,  chemical,  and  biological  mechanisms.  For  example,  dispersion  can  increase 
the  effectiveness  of  adsorption  because  "spreading"  of  the  contaminant  plume  exposes  the 
pollutants  to  a  greater  volume  of  soil  particles  at  a  reduced  concentration.  In  addition, 
dispersion  increases  the  potential  for  biological  reactions  which  rely  upon  exposure  to 
nutrients,  bacteria  and/or  oxygen.  Information  about  chemical  mechanisms,  although 
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important,  has  not  been  particularly  useful  to  fate  and  transit  modelers.  Biodegradation  is 
currently  receiving  the  most  attention  as  a  remediation  option,  primarily  due  to  the  advantages 
associated  with  this  technology. 

Although  natural  attenuation  is  an  attractive  option,  it  will  not  be  appropriate  for  every 
site.  One  must  consider  many  factors  when  choosing  a  remediation  strategy.  As  much 
information  ^ould  be  gathered  as  possible  in  supporting  a  decision. 

Models  have  been  designed  to  assist  with  predictions  about  the  fate  and  transport  of 
contaminants.  BIOPLUME  II™  is  one  of  those  models  and  it  is  used  in  this  study  to  gain  a 
better  understanding  of  the  parameters  and  mechanisms  affecting  natural  attenuation  in  the 
saturated  2one. 
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III. 


Introduction 

From  the  review  of  the  numerous  mechanisms  and  parameters  associated  with  the 
natural  attenuation  process  in  the  last  chapter,  it  is  evident  that  the  development  of  a  model 
which  includes  every  mechanism  and  parameter  would  be  a  formidable  task  indeed.  In  fact, 
no  such  model  exists. 

When  using  models  to  predict  the  fate  and  transport  of  a  contaminant  plume,  one  must 
be  conscious  of  the  limitations  of  the  model  and  in  light  of  those  limitations,  view  the  output 
with  caution.  In  this  study,  the  objective  is  to  examine  some  of  the  parameters  which  affect 
natural  attenuation  and  to  formulate  some  general  guidelines  which  will  indicate  the  feasibility 
of  natural  attenuation  as  a  remediation  alternative  for  given  conditions. 

Reaeach  pgaga 

To  meet  this  research  objective,  the  computer  model,  BIOPLUME  II™  will  be 
employed  to  simulate  the  fate  and  transport  of  a  theoretical  contaminant  plume  as  it  migrates 
through  a  theoretical  site  under  various  site  conditions. 

Rifai  and  others  have  conducted  a  sensitivity  analysis  of  the  model.  This  sensitivity 
analysis  has  been  included  as  Appendix  E  for  the  convenience  of  the  field  engineer 
investigating  the  effects  of  these  parameters  on  natural  attenuation.  From  the  sensitivity 
analysis,  it  was  concluded  that  the  biodegraded  mass  in  the  model  is  most  sensitive  to 
hydraulic  conductivity  and  the  coefficients  of  reaeration  and  anaerobic  biodegradation.  Tliis 
sensitivity  analysis  focused  primarily  on  the  percentage  of  mass  biodegraded  and  ignored  the 
distance  that  the  contaminant  plume  migrates.  Depending  on  the  characteristics  of  a  given 
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site,  the  distance  a  contaminant  plume  is  expected  to  migrate  may  be  more  relevant  than  the 
amount  of  contaminant  mass  biodegraded.  The  current  study  consists  of  two  parts.  The  first 
part  will  repeat  the  sensitivity  analysis  conducted  by  Rifai  with  some  minor  modifications  and 
focus  on  the  distance  which  the  contaminant  plume  migrates  rather  than  the  amount  of 
contaminant  mass  biooegraded.  Additionally,  this  study  will  further  examine  the  three 
parameters  identified  as  being  the  most  significant  to  biodegradation  by  varying  all  three  over 
a  range  of  values  rather  than  one  at  a  time,  as  is  typically  done  in  a  sensitivity  analysis.  By 
varying  all  three  parameters,  a  general  relationship  can  be  postulated  regarding  the  effect  of 
these  parameters  on  natural  attenuation. 

Pan  I  -  Sensitivity  Analysis 

BIOPLUME  n™  will  be  applied  to  theoretical  plumes  of  JP-4  of  constant  source 
strength  and  duration.  The  model  will  be  set  up  to  analyze  the  natural  attenuation  of  the 
contaminant  plume  and  to  project  the  distance  the  plume  will  travel  in  a  specified  length  of 
time. 

To  estimate  which  parameters  have  the  most  significant  effect  on  the  distance  a 
contaminant  plume  travels,  a  sensitivity  analysis  will  be  conducted.  Natural  attenuation 
mechanisms  included  in  the  model  are  dispersion,  sorption,  and  biodegradation.  Both  aerobic 
and  anaerobic  biodegradation  can  be  simulated  (28:  2-5). 

To  begin  the  analysis  of  the  parameters  affecting  natural  attenuation  in  the  saturated 
zone,  a  hypothetical  site  will  be  simulated.  The  hypothetical  site  will  be  a  simple, 
homogeneous  aquifer  with  constant  hydraulic  gradient  and  a  single  contaminant  source.  The 
site  is  mapped  as  an  eleven-cell  by  twenty-cell  grid.  The  initial  site  dimensions  will  be  500 
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feet  by  1000  feet.  A  simulated  injection  well  will  serve  as  the  contaminant  source  and  five 
simulated  observation  wells  will  be  strategically  located  downgradient  from  tte  injection  well 
at  regular  intervals  to  track  the  migration  of  the  contaminant  plume  through  the  site.  The  ate 
and  wells  are  depicted  in  Figure  6.  Conattering  diat  for  the  majority  of  reaoration  projects, 
the  contaminant  source  is  removed,  the  injection  well  simulation  should  serve  as  a  wora  case 
scenario  and  yield  conservative  guidance. 

The  model  output  shows  a  maximum  concentration  for  the  cell  in  which  the  source  is 
located.  This  concentration  is  a  maximum  average  concentration  in  the  SO  ft  by  SO  ft  cell  and 
not  the  concentration  precisely  at  the  source.  The  diaance  is  also  not  precise  but  related  to 
the  cell  size,  in  this  case  SO  ft.  The  diaance  that  a  plume  travels  in  SO  years  is  considered  to 
be  the  incremental  diaance  to  the  fira  unaffected  grid  cell. 

Input  variables  for  the  BIOPLUME  II™  model  are  liaed  in  Appendix  D.  Values  for 
the  variables  of  interea  are  obtained  from  the  literature.  Table  4  lias  initial  site  data  which  is 
pertinent  to  this  audy  and  was  used  for  the  base  run. 

In  order  to  make  comparisons  with  Rifai's  work,  the  same  parameters  were  chosen  to 
be  inveaigated.  These  parameters  include,  porosity,  hydraulic  conductivity,  and  longitudinal 
dispersivity.  Because  JP-4  is  less  dense  than  water  and  floats  on  top  of  the  water  table, 
reaeration  by  oxygen  transfer  from  the  vadose  zone  may  be  potentially  significant  to  natural 
attenuation  of  JP  -  4  (38:4).  Rifai  also  found  these  parameters  to  be  significant  contributors 
to  biodegradation.  Therefore,  a  reaeration  parameter  will  be  introduced  as  well  as  a  rate  of 
decay  term  to  simulate  anaerobic  biodegradation. 
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Figure  6.  Simulated  Site  and  Observation  Weils 


Table  4.  Input  data  for  initial  computer  simulation 


Simulation  time 

SO  years 

Grid  Size 

11  X  20 

Cell  Size 

so  ft  X  50  ft 

Poroaty 

3 

Longitudinal  Dispersivity 

10  ft 

Transverse  Dispersivity 

3  ft 

Txx 

.0025  ft*/s 

Tyy 

.0025  ft*/s 

Aquifer  Thickness 

25  ft 

Hydraulic  Gradient 

4.29E-3  fl/ft 

Injection  Well  at  Cell 

maamam 

Injection  Rate 

0.0002  cfs 

Cone,  of  Contaminant  in  injected  water 

ISO  mg/1 

Cone,  of  Oxygen  in  Injected  Water 

0.0  mg/1 

Cone,  of  Natural  Recharge  of  Oxygen 

8.0  mg/1 
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Porosity 

Typical  values  for  porosity  were  listed  in  chapter  two.  Porosity  values  will  be  varied 
from  .1  to  .6.  The  variation  in  the  distance  a  contaminant  plume  travels  over  this  range  of 
porosities  will  be  shown.  Since  velocity  is  inversely  proportional  to  porosity,  it  is  expected 
that  the  distance  a  plume  travels  in  a  given  time  period  will  be  significantly  afiected. 
Hydraulic  Conductivity 

Some  values  for  hydraulic  conductivity  were  also  given  in  chapter  two.  There  are 
some  model  limitations  which  prevent  examining  the  full  range  of  hydraulic  conductivities 
presented  in  chapter  two;  however,  a  trend  should  be  evident  from  the  examined  range  from 
which  general  guidelines  can  be  formulated.  Because  there  is  no  direct  input  parameter  in  the 
model  for  hydraulic  conductivity,  the  transmissivity  and  aquifer  thickness  will  be  varied  to 
vary  the  hydraulic  conductivity  from  10'^  to  10"*  fi/sec.  Hydraulic  conductivity  is  related  to 
transmissivity  and  thickness  by  the  following  relationship; 


T 

K=  — 

b 


where 


K  -  hydraulic  conductivity  (ft/sec) 

T  -  transmissivity  (ftVsec) 
b  -  aquifer  thickness  (ft) 

The  distance  the  contaminant  plume  travels  will  be  plotted  over  this  range  of  values.  Rifai's 
analysis  showed  hydraulic  conductivity  to  have  a  significant  efiect  on  biodegradation  and 
from  the  above  relationship  it  is,  likewise,  expected  to  affect  plume  transport. 
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Dispersivity 

Values  for  longitudinal  dispersivity  will  be  varied  from  1  ft  to  10  ft  and  the 
relationship  to  plume  travel  plotted.  Because  longitudinal  dispersion  spreads  the  contaminant 
in  the  direction  of  groundwater  flow,  it  is  expected  that  longitudinal  dispersivity  will  have  an 
effect  on  the  distance  the  plume  travels  in  50  years. 

Reaeration 

The  reaeradon  coefficient  will  be  varied  over  the  same  range  which  Rifai  used,  0.001 
per  day  to  0.005  per  day.  A  plot  will  present  the  results.  From  Rifai's  sensitivity  analysis,  it 
is  expected  that  the  reaeration  coefficient  will  have  a  significant  effect  on  the  distance  a 
contaminant  plume  travels. 

Anaerobic  Degradation 

Anaerobic  degradation  wUl  be  evaluated  by  varying  the  rate  of  decay  from  0.002  per 
day  to  .01  per  day.  According  to  Rifai,  anaerobic  degradation  can  be  expected  to  have  a 
significant  effect  on  the  contaminant  plume. 

Results  from  the  simulations  will  be  compiled  and  presented  in  grt^hical  format 
showing  the  effect  of  varying  different  parameters  on  maximum  concentration  and  plume 
travel.  The  output  data  from  the  model  should  show  how  the  value  of  various  parameters 
affect  the  movement  of  the  contaminant  plume,  contaminant  levels  along  the  plume,  and  the 
time  and  distance  that  it  takes  the  contaminant  to  degrade  to  acceptable  limits.  These 
relationships  will  be  illustrated  by  plotting  them  on  a  graph  which  will  show  how  changes  in 
the  value  of  a  parameter  affect  these  plume  characteristics. 
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Pan  n 


The  ate  and  site  conditions  to  be  simulated  in  the  second  pan  of  tin  audy  were  the 
same  as  the  base  conditions  for  Pan  1. 

The  reaeradon  coefficient  was  varied  over  a  range  from  .001  day  ’  U)  .005  day  '  and 
the  coefficient  of  anaerobic  decay  varied  from  .002  day  '  to  .01  day  '. 

For  the  first  set  of  data,  the  hydraulic  conductivity  was  held  constant  at  lO**  fr/sec  as 
the  values  for  reaeration  and  anaerobic  decay  were  varied  over  their  respective  ranges.  This 
approach  produced  25  sets  of  data  input  for  the  model  for  each  value  of  hydraulic 
conductivity.  This  procedure  was  repeated  for  hydraulic  conductivity  values  of  10'^  ft/sec, 

10*^  ft/sec,  and  finally,  for  10'^  ft/sec  for  a  total  of  100  data  sets.  Computer  simulations  for 
hydraulic  conductivity  of  10'^  were  rather  time  consuming  and  values  greater  than  10'^  ft/sec 
would  have  been  infeasible  to  investigate  in  this  study.  A  regression  analyas  was  performed 
on  this  data  to  generate  a  mathematical  relationship  between  these  parameters. 

Smnthjffy 

The  research  design  and  analysis  methodology  that  were  used  in  this  research  project 
were  designed  to  examine  several  of  the  parameters  involved  in  the  natural  attenuation  of  JP-4 
in  the  saturated  zone.  The  effects  that  the  three  parameters:  hydraulic  conductivity, 
reaeration,  and  anaerobic  biodegradation  have  on  the  natural  attenuation  of  JP-4  were  further 
studied  to  estimate  how  these  three  parameters  in  conjunction  with  one  another  affect  natural 
attenuation. 
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IV. 


Introduction 

The  analyas  of  the  results  from  this  research  are  presented  in  the  order  outlined  in 
chapter  three.  Part  one,  a  sensitivity  analysis  of  the  parameters  pmoaty,  hydraulic 
conductivity,  longitudinal  diqiersivity,  reaenuion,  and  anaerobic  biodegradation  will  be 
presented  followed  by  a  discussion  of  pan  two,  the  sensitivity  and  regression  analysis  of 
hydraulic  conductivity,  reaeration,  and  anaerobic  biodegradation  varying  in  concert 
Pan  1  -  Senativity  Analyaa 

Porosity 

In  order  to  investigate  the  effect  porosity  has  on  natural  attenuation,  the  values  were 
varied  as  described  in  charter  three  while  holding  all  other  baseline  parameters  constant 
Porosity  was  varied  from  .1  to  .6,  covering  the  full  range  of  porosities  presented  in  chapter 
two.  The  results  are  presented  in  Table  5  and  Figure  7.  Rifai  concluded  from  her  sensitivity 
analysis  that  porosity  did  not  have 

a  agnificant  impact  on  biodegradation;  however,  these  results  indicate  that  porosity  does  have 
a  significant  impact  on  the  distance  which  a  contaminant  plume  is  expected  to  travel  and  an 
unpredictable  effect  on  maximum  concentration  levels.  Combining  the  results  of  the  two 
analyses,  one  may  be  able  to  predict  that  increasing  porosity  decreases  the  distance  the  plume 
will  migrate  in  a  given  time  period,  but  porosity  does  not  have  a  predictable  effect  on 
biodegradation. 
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Table  5  Data  for  varying  porosity 


Well  # 

Hydrocarbon 

OxvQen  1 

Max  Cone. 

Distance 

time  to  reach  well  (years) 

mq/liter 

Feet 

Source 

0 

0 

53 

0 

1 

1.633 

1  02 

46.2 

250 

2 

2  449 

1  429 

43.2 

300 

3 

3.061 

2  041 

40.2 

350 

4 

3.878 

2.653 

37.9 

400 

5 

4694 

3469 

35  8 

450 

Porositv  .2 

Source 

0 

0 

55 

0 

1 

3.252 

2.033 

46 

250 

2 

4  878 

2.846 

42.7 

300 

3 

6.098 

4.065 

40 

350 

4 

7.724 

5.285 

37.7 

400 

5 

9.35 

6.911 

357 

450 

Porositv  .3 

Source 

0 

0 

53 

0 

1 

4.878 

3  049 

46 

250 

2 

7.317 

4.26S 

42.7 

300 

3 

9.146 

6.093 

40 

350 

4 

11.585 

7.927 

37.7 

400 

5 

14.024 

10.366 

35.7 

450 

Poros 


Source 

0 

0 

56 

0 

1 

6.452 

4.032 

45.6 

250 

2 

9.677 

5.645 

42.9 

300 

3 

12.903 

8.085 

40.5 

350 

4 

15.323 

10.484 

38.2 

400 

5 

18.548 

13.71 

35.8 

450 

Porositv  .5 

Source 

0 

0 

61 

0 

1 

8.163 

5.102 

45.9 

250 

2 

12.245 

7.143 

43 

300 

3 

15.306 

10.204 

40.1 

350 

4 

19.388 

13.265 

37.8 

400 

5 

23.469 

17.347 

35.7 

450 

F>orositv  .6 

Source 

0 

0 

55 

0 

1 

9.756 

6.098 

45.5 

250 

2 

14.634 

8.537 

42.7 

300 

3 

18.293 

12.195 

39.8 

350 

.  4 

23.171 

15.854 

37.4 

400 

5 

28.049 

23.171 

34.6 

450 
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Distance  (Feet)  Maximum  Concentration  (mg/L) 


Porosity 
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Figure  7.  Changes  in  maximum  concentratimi  and  distance  widi  changes  in  pofosiqr 
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Hydraulic  Conductivity 

Hydraulic  Coi^uctivity  was  varied  from  10*^  to  10'^  fr/sec.  The  results  are  {nesent^ 
in  Table  6  and  Figure  8  and  Figure  9.  Hydraulic  Conductivity  had  the  most  significant  efrect 
of  any  parameter  on  the  behavior  of  the  contaminant  plume.  Both  maximum  contaminant 
concentrations  and  plume  migration  distance  were  agnificantly  affected  over  the  range  of 
hydraulic  conductivity  values.  Because  of  the  sgnifrcant  impact  this  parameter  has  on  plume 
behavior,  it  was  chosen  as  one  of  the  three  parameters  to  be  investigated  in  Part  n  of  this 
study.  The  results  of  the  sensitivity  analysis  show  that  the  larger  the  value  of  the  hydraulic 
conductivity  at  a  site,  the  less  favorable  conditions  are  for  natural  attenuation. 

Longitudinal  Dispersivity 

Longitudinal  EMspersivity  was  varied  from  1  ft.  to  10  ft  The  results  are  presented  in 
Table  7  and  Figure  10.  Longitudinal  Dispersivity  in  this  range  had  negligible  effect  on  the 
distance  a  contaminant  plume  is  expected  to  migrate.  There  was  a  slight  effect  on  the 
maximum  contaminant  concentration.  From  both  the  percent  mass  biodegraded  in  Rifai's 
sensitivity  analysis  and  the  distance  which  a  contaminant  plume  is  expected  to  migrate,  it  can 
be  concluded  that  longitudinal  disperavity  does  not  have  a  significant  impact  on  natural 
attenuation  of  JP-4  in  the  saturated  zone. 

Reaeration 

The  reaeration  coefficient  was  varied  from  .001  day  ’  to  .005  day‘‘  .  The  results  are 
presented  in  Table  8  and  Figure  11.  It  is  evident  from  the  results  that  reaeradon  has  a 
significant  effect  on  both  the  diaance  which  a  contaminant  plume  is  expected  to  migrate  and 
on  the  maximum  contaminant  concentration  in  the  plume.  Rifai  also  concluded  that 
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Table  6.  Data  for  varying  hydraulic  conductivity 


drauiic 


Well  # 


Source 

1 

2 

3 

4 

5 


Hydraulic 


Source 

1 

2 

3 

4 

5 


Hydraulic 


Conductivity  10E-06 


Hydrocarbon  O 


time  to  reach  well  (years 


0  0 
>50  >50 


Max  Cone  Distance 


feet 


131  0 

0  250 

300 
350 
400 
450 


Conductivity  10E-05 


0  0 

27.273  11.364 

43.1 02  25 

>50  40.909 

>50 


Conductivity  10E-04 


150 

90.6 

6.6 
0 


0 

250 

300 

350 

400 

450 


I  Source 

0 

0 

53 

0 

j  1 

4.878 

3.049 

46 

250 

2 

7.317 

4.268 

42.7 

300 

3 

9.146 

6.098 

40 

350 

4 

11.585 

7.927 

37.7 

400 

5 

14.024 

10.366 

35.7 

450 

!  Hydraulic  Conductivity  10E-03 

Source 

0 

0 

9 

0 

1 

31.579 

21.053 

0.7 

5000 

2 

44.737 

31.579 

0.2 

7000 

>50 

42.105 

0 

9000 

50 

11000 

1 

i  5 

>50 

13000 

Maximum  Concentration  (mg/l) 


Hydraulic  Conductivity 


0  ^ ^ ^ ^ ^ - — 

1E-06  1E-05  0.0001  0.001 

Hydraulic  Conductivity  (ft/sec) 

1 


Figure  8.  Changes  in  maximum  concentration  wiA  changes  in  hydraulic  conductivity 
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Distance  (feet) 


Hydraulic  Conductivity 


12000 


10000  - 


8000  f 


6000 


4000 


2000 


IE-06  IE-05  0.0001  0.001 

Hydraulic  Conductivity  (ft/sec) 


Figure  9.  Changes  in  distance  widi  changes  in  hydraulic  conductiviQr 
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Table  7  Data  for  varying  longitudinal  dispersivity 


Loqitudinal  Dispersivity  1  ! 

Well  # 

Hydrocarbon 

i  Oxvoen 

Max  Cone.  I 

Distance 

time  to  reach  well  (years) 

■Mian 

Feet 

Source 

0 

0 

57 

0 

1 

6.093 

4.235 

56.4 

250 

2 

9.146 

6.098 

58.7 

300 

3 

10.976 

6.537 

58.4 

350 

4 

13.415 

10.976 

57.8 

400 

5 

15.854 

13.415 

56.7 

450 

Loqitudinal  Dispersivitv  2 

Source 

0 

0 

56 

0 

1 

6.098 

3.659 

569 

250 

n 

4- 

7.927 

6.098 

56.5 

300 

3 

10.366 

7.927 

55.2 

350 

4 

12.805 

10.366 

54.2 

400 

5 

15  854 

12.805 

53.2 

450 

ituainal  Dis 

ersiv 

Source 

1 

2 

3 

4 

5 


3 


0 

5  488 

7.927 

10.366 

12.805 

15.244 


itudina!  Dispersivitv  4 


0 

3.659 

5.488 

7.927 

9.756 

12.195 


56 

55.2 

54.3 

52.3 

51.3 
5 


Source 

1 

2 

3 

4 

5 

0 

5.488 

7.927 

10.366 

12.805 

15.244 

0 

3.659 

5.488 

7.927 

9.756 

12.195 

56 

53.6 
52.3 

49  8 

48.6 

47.1 

0 

250 

300 

350 

400 

450 

Loqitudinal  Dispersivitv  5 

I  Source 

0 

0 

55 

0 

^  1 

5.488 

3.049 

52 

250 

I  2 

7.317 

4.878 

50.4 

300 

3 

9.756 

7.317 

47  8 

350 

4 

12.195 

9.146 

46.2 

400 

5 

14.634 

11.585 

44.6 

450 

Table  7  (cent.)  Data  for  varying  longitudinal  dispersivity 


■ILoqitudinai  Dispersivitv  6 


.  Source 

0 

55 

0 

1 

4.878 

3.049 

50.6 

250 

'  2 

7.317 

4.878 

48.6 

300 

3 

9.756 

6.707 

45.9 

^ 

1  1 
Ww  W 

4 

12.195 

9.146 

44 

430 

5 

14.634 

10.976 

42.6 

450 

i*Loaitudinal  0 

ispersivitv  7 

Source 

0 

0 

54 

0 

I  1 

4.878 

3.049 

49.4 

250 

2 

7.317 

4.876 

^6  5 

300 

'  3 

9.756 

6.707 

44.3 

350 

;  4 

12.195 

8.537 

42.2 

400 

5 

14.024 

10.976 

40.5 

450 

ilLooitudinal  Dispersivitv  8  li 

1  Source 

0 

0 

54 

0  ! 

i  1 

4.878 

3.049 

48.2 

250  I 

t  2 

7.317 

4.876 

45.4 

3 

9.146 

6.707 

42.8 

350  1 

4 

12.195 

8.537 

40.6 

400  j 

5 

14.024 

10.366 

38.9 

450  1 

Loaitudinal  D 

ispersivitv  9  II 

Source 

0 

0 

53 

0 

1 

4.878 

3.049 

47.1 

250 

2 

7.317 

4.268 

44 

300 

3 

9.146 

6.098 

41.4 

350 

4 

12.195 

8.537 

39.1 

400 

5 

14.024 

10.366 

37.2 

450  1 

Loaitudinal  D 

ispersivitv  10  I 

Source 

0 

0 

53 

0 

!  1 

4.878 

3.049 

46 

250 

2 

7.317 

4.258 

42.7 

300 

3 

I 

9.146 

6.098 

40 

350 

4 

11.585 

7.927 

37.7 

400 

5 

14  024 

10.366 

35  7 

450  i 

46 


Distance  (Feet)  Maximum  Concentration  (mg/I) 


Longitudinal  Dispersivity 


Long.  Dispersivity  (Feet) 


Long.  Dispersivity  (Feet) 


Figure  10.  Changes  in  maximum  concentration  and  distance  with  changes  in  longitudinal 
diq)eTsivity 
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Table  8.  Data  for  varying  reaeration  coefficient 


jReaeration  Coeff 

cient  uu  1 

— ^ ^ - 

1 

il  Wei!  # 

I  Oxvaen 

!  Max  Cone. 

i  Distance 

- 

ii 

time  to  reach  well  (years) 

mg/liter 

feet 

Source 

0 

0 

43 

0 

1 

5.4S8 

3.049 

13 

250 

2 

7.927 

4.266 

7.6 

300 

3 

11.585 

6.707 

4.3 

350 

! 

4 

-1 7. 073 

8  537 

2,4 

400 

1 

1 

_ 5 _ 

2S.659 

10.976 

1.3 

450 

1 

1 

-  i 

ilReaeration  Coefficient  .002  ll 

Source 

0 

0 

39 

0 

ii 

it 

1 

6.098 

3.049 

4.2 

OCA 

•• 

2 

26  829 

4.876 

1.5 

300 

I' 

3 

-“50 

6.707 

0.3 

350 

i. 

4 

9  14'.' 

0 

400 

ii 

i  K 

I  ^ 

-trt  Qyc> 

i  w.  9  r  w 

0 

450 

iReaeratior  Coeffi 

w  t  w  •  4  w  w 

"ii 

t, 

J; 

1  oource 

0 

0 

35 

r 

w 

i> 

;  1 

26.820 

3.049 

1,4 

250 

i< 

2 

>50 

4.S78 

0.1 

300 

f 

i  3 

7.317 

0 

350 

4 

9.146 

0 

400 

i  5 

11.585 

0 

450 

‘iReaeration  Coefficient  .004  |! 

Source 

0 

0 

33 

0 

1 

1 

28.659 

3.049 

0.4 

250 

2 

>50 

4.878 

0 

300 

1 

3 

7.317 

0 

350 

1 

4 

9.146 

0 

400 

5 

11.585 

0 

450 

jReaeration  Coefficient  .005  II 

j  Source 

0 

0 

31 

0 

1- 

>50 

3.049 

0 

250 

■ 

2 

4.878 

0 

300 

1 

!  3 

7.317 

0 

350 

‘  4 

9146 

0 

400 

■ 

!  5 

11.585 

0 

450 

i 
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Reaeration 
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0.001  0.002  0.003  0.004  0.005 

Reaeration  Coefficient  (l/day) 


Reaeration  Coeffident  (l/dayH 

I 


Figure  1 1.  Granges  in  maximum  concentration  and  distance  with  changes  in  reaeration 
coefficient 
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leaeration  had  a  significant  effect  on  biodegradation  of  a  contaminant  plume.  Because 
reaeration  appears  to  be  a  significant  contributor  to  the  natural  attenuation  process,  it  was 
included  as  a  parameter  of  interest  in  Part  n  of  this  study.  Results  from  ttese  analyses 
indicate  that  increasing  reaeration  is  important  to  natural  attenuation  and  the  higher  the 
coefficient  of  reaeration  the  mote  likely  it  is  that  natural  attenuation  will  be  a  candidate  for  a 
remediation  of  the  site. 

Anaerobic  Biodegradation 

Values  for  the  coefficient  of  anaerobic  biodegradation  were  varied  from  .002  day  '  to 
.01  day  '.  The  results  are  presented  in  Table  9  and  Figure  12.  It  was  anticipated  that 
anaerobic  biodegradation  would  have  a  similar  effect  as  reaeration  because  both  are  handled 
by  the  model  in  a  similar  manner.  This  was  indeed  the  case;  however,  anaerobic 
biodegradation  appears  to  have  a  slightly  more  significant  effect  than  reaeration.  This  result 
is  probably  due  to  the  fact  that  the  contaminant  concentration  near  the  source  or  near  the 
center  of  the  plume  increases  as  the  oxygen  at  the  center  of  the  plume  is  being  depleted  by 
the  aerobic  biodegradation  process.  Once  all  oxygen  is  completely  gone  and  anaerobic 
microbes  become  acclimated,  anaerobic  biodegradation  progresses  faster  than  aerobic 
biodegradation.  As  was  noted  in  chapter  two,  however,  anaerobic  biodegradation  is  less 
certain  and  the  anaerobic  microbes  take  time  to  acclimate.  The  model  does  not  account  for 
the  acclimation  period. 

Wilson  has  measured  anaerobic  biodegradation  rates  of  .002  day  '  for  certain  chemicals 
at  one  site  (37:8).  The  model  indicates  that  even  rates  this  low  may  have 
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Table  9.  Data  for  varying  coefficient  of  anaerobic  biodegradation 


Anaerobic  Coefficient  .002  i 

Well# 

Hydrocarbon 

Oxvaen  i 

Max  Cone. 

Distance 

time  to  reach  well  (years) 

mq./iiter 

Feet 

Source 

n 

w 

0 

39 

0 

1 

6.098 

3.049 

4.2 

250 

2 

26.829 

4.858 

1.5 

300 

3 

31.707 

7.317 

0.3 

350 

4 

>50 

9.146 

0 

400 

5 

10.976 

0 

450 

Anaerobic  Coeffi 

cient.004 

ilAnaerobic  Coefficient.008 


0 

049 
878 
7.317 
9.146 
11.585 


Anaerobic  Coefficient  .006 

Source 

0 

0 

30 

0 

1 

>50 

3.049 

0 

250 

2 

5.488 

0 

300 

3 

7.927 

0 

350 

4 

9.756 

0 

400 

5 

12.195 

0 

450  i 

Source 

0 

0 

28 

0 

1 

>50 

3.049 

0 

250 

2 

5.488 

0 

300 

3 

7.927 

0 

350 

4 

9.756 

0 

400 

5 

12.195 

0 

450 

Anaerobic  Coefficient  .01  H 

Source 


Anaerobic  Degradation 


Coefficient  of  Anaerobic  Decay  (1/day) 


Anaerobic  Degradation 


Coefficient  of  Anaerobic  Decay  (l/day) 

Figure  12.  Changes  in  maximum  concentration  and  distance  with  changes  in  coefficient  of 
anaerobic  biodegradation 
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significant  impact  on  natural  attenuation.  If  anaerobic  degradation  can  be  verified  at  a  ate,  it 
can  be  concluded  that  it  has  a  significant  effect  on  natural  attenuation. 

EadJL 

The  three  parameters  chosen  for  investigation  were  hydraulic  conductivity,  coefficient 
of  reaeration,  and  anaerobic  degradation.  These  three  parameters  were  chosen  because  they 
were  identified  in  both  Rifai's  sensitivity  analysis  and  the  Part  I  senativity  analysis  as  having 
a  significant  impact  on  contaminant  plume  behavior. 

The  Part  II  study  was  conducted  to  derive  a  general  relationship  between  these 
parameters  as  they  relate  to  natural  attenuation.  The  results  from  100  computer  simulations 
are  shown  in  tables  10  -  13.  The  tables  riiow  the  value  for  each  parameter,  the  distance 
which  a  contaminant  plume  is  expected  to  migrate  in  SO  years  and  the  maximum 
concentration  at  the  source  for  the  various  combinations  of  parameter  values.  As  can  be  seen, 
hydraulic  conductivity  was  the  dominant  parameter.  Neither  distance  nor  maximum 
contaminant  concentrations  were  significantly  afiected  by  variations  in  the  coefficient  of 
leaeration  and  anaerobic  decay  for  a  given  hydraulic  conductivity.  For  a  hydraulic 
conductivity  value  of  10'^  ft/sec,  expected  migration  distance  of  the  contaminant  plume 
begins  to  be  affected  when  the  other  two  parameters  are  varied. 
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Table  10:  Computer  Simulations  for  Hydraulic  Conductivity  of  10'*  ft/sec 


RC 

ADC 

MAX  CONC 

DISTANCE 

.001 

.002 

3 

150  -  200 

.001 

.004 

3 

150  -  200 

.001 

.006 

3 

150  -  200 

.001 

.008 

3 

150  -  200 

.001 

.010 

3 

250 

.002 

.002 

3 

150 

.002 

.004 

3 

150 

.002 

.006 

3 

150 

.002 

.008 

3 

150 

.002 

.010 

3 

150 

.003 

.002 

3 

150 

.003 

.004 

3 

150 

.003 

.006 

3 

150 

.003 

.008 

3 

150 

.003 

.010 

3 

150 

.004 

.002 

3 

150 

.004 

.004 

3 

150 

.004 

.006 

3 

150 

.004 

.008 

3 

150 

.004 

.010 

3 

150 

.005 

.002 

3 

150 

.005 

.004 

3 

150 

.005 

.006 

3 

150 

.005 

.008 

3 

150 

.005 

.010 

3 

150 
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Table  11:  Computer  Simulations  for  Hydraulic  Conductivity  of  10'^  ft/sec 


RC 

ADC 

MAX  CONC 

DISTANCE 

.001 

.002 

20 

250 

.001 

.004 

19 

250 

.001 

.006 

19 

250 

.001 

.008 

19 

250 

.001 

.010 

19 

250 

.002 

.002 

19 

250 

.002 

.004 

19 

250 

.002 

.006 

19 

250 

.002 

.008 

19 

250 

.002 

.010 

19 

250 

.003 

.002 

19 

250 

.003 

.004 

19 

250 

.003 

.006 

19 

250 

.003 

.008 

19 

250 

.003 

.010 

19 

250 

.004 

.002 

19 

250 

.004 

.004 

19 

250 

.004 

.006 

19 

250 

.004 

.008 

19 

250 

.004 

.010 

19 

250 

.005 

.002 

19 

250 

.005 

.004 

19 

250 

.005 

.006 

19 

250 

.005 

.008 

19 

250 

.005 

.010 

19 

250 
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Table  12.  Computer  Simulations  for  Hydraulic  Conductivity  of  lO'*  ft/sec 


RC 

ADC 

MAX  CONC 

DISTANCE 

.001 

.002 

35 

1800 

.001 

.004 

31 

2000 

.001 

.006 

29 

2000 

.001 

.008 

27 

2000 

.001 

.010 

26 

2000 

.002 

.002 

33 

2000 

.002 

.004 

30 

2000 

.002 

.006 

28 

2000 

.002 

.008 

27 

2000 

.002 

.010 

26 

2000 

.003 

.002 

31 

2000 

.003 

.004 

29 

2000 

.003 

.006 

27 

2000 

.003 

.008 

26 

2000 

.003 

.010 

25 

2000 

.004 

.002 

30 

2000 

.004 

.004 

28 

2000 

.004 

.006 

27 

2000 

.004 

.008 

26 

2000 

.004 

.010 

25 

2000 

.005 

.002 

29 

2000 

.005 

.004 

27 

2000 

.005 

.006 

26 

2000 

.005 

.008 

25 

2000 

.005 

.010 

25 

2000 
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Table  13;  Computer  Simulations  for  Hydraulic  Conductivity  of  10^  ft/sec 


Rsgreasfln 

To  derive  the  relationship  between  these  parameters,  as  they  relate  to  the  natural 
attenuation  process  of  JP-4  in  the  saturated  zone,  a  regression  analysis  was  performed  on  die 
data  from  these  computer  simulations.  The  regression  analysis  suggested  the  following 
mathematical  relationships; 


S  =  880.2644  +  2133161.421(HC)  -  21750  (RC)  -  20750  (ADC) 
C  =  18.6688  -10994.81377  (HC)  +  210  (RC)  -  232.5  (ADC) 


where 


S  -  Instance,  in  feet,  which  the  plume  is  expected  to  travel  in  50  years 
C  -  Maximum  average  concentration,  in  mgAiter,  in  the  cell  in  which  the 
injection  well  is  located 
HC  -  Hydraulic  Conductivity  (ft/sec) 

RC  -  Reaeration  Coefficient  (day  *) 

ADC  -  Anaerobic  Decay  Coefficient  (day  ') 


It  is  evident  from  the  parameter  coefficients  in  these  equations  how  signiricant  hydraulic 
conductivity  is  to  plume  behavior.  Statistical  results  indicate  that  the  distance  equation  is 
more  useful  than  the  concentration  equation.  The  R^  value  for  the  distance  equation  was 
.6448  and  the  R^  value  for  the  concentration  equation  was  .1981.  The  R^  measures  the 
goodness  of  rit  of  a  regression  model.  It  measures  the  proportion  of  variance  of  the 
regression  analysis. 

The  field  engineer  may  be  able  to  use  these  relation^ips  to  gain  a  precursory  idea 
whether  or  not  natural  attenuation  should  be  considered  as  a  remediation  option.  If  the 
distance  to  an  exposure  point  is  known,  the  distance  equation  can  be  used  to  estimate  how  far 
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the  plume  is  expected  to  travel  in  SO  years.  If  the  calculated  distance  is  much  less  than  the 
distance  to  the  exposure  point,  natur?l  attenuation  may  suffice  and  more  costly  remediation 
technologies  may  not  need  to  be  implemented.  The  concentration  equation  can  be  used  by 
calculating  the  concentration  at  the  source.  If  the  calculated  concentration  is  much  less  than 
regulatory  limits,  further  remediation  may  not  be  necessary. 

Some  General  Observations 

Tile  BIOPLUME  II™  model  exhibited  some  unusual  characteristics.  Hiese 
characteristics  may  have  been  due  to  simplifications  made  in  this  study  or  may  be  inherent  to 
the  model.  The  model  seemed  to  generate  a  "pulang"  contaminant  plume  when  a  reaeration 
coefficient  or  coefficient  of  anaerobic  decay  are  introduced. 

An  article  by  David  L.  Freyburg  in  a  1988  issue  of  Groundwater  magazine  suggests 
that  calibration  of  a  model  does  not  necessarily  lead  to  accurate  predictions.  Also,  a  model  is 
only  as  good  as  the  data  you  put  into  it.  In  light  of  these  cautions,  it  would  be  unwise  to 
extend  the  results  of  this  study  alone  to  a  site  and  decide  to  implement  natural  attenuation 
based  on  these  findings.  Hopefully,  these  results  have  provided  some  general  guidelines 
which  can  be  followed  in  choosing  a  remediation  alternative. 
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Introduction 

In  reviewing  the  literature  and  previous  studies  of  natural  attenuation  mechanisms,  it 
was  evident  that  the  study  of  natural  attenuation  is  new  and  not  yet  well  understood.  There  is 
extensive  ongoing  research  to  gain  a  better  understanding  of  the  mechanisms  and  the 
parameters  involved  in  the  natural  attenuation  process.  As  evidenced  in  this  study,  the 
number  of  mechanisms  and  parameters  are  numerous  indeed.  Coupled  with  the  number  of 
mechanisms  and  parameters  are  the  differing  site  conditions  and  numerous  chemicals  of 
concern. 

MamraLAtte-nuation  of  JP-4  in  ibc-Satigatsd.?Qnc 

It  is  well  recognized  that  certain  chemical  contaminants,  including  JP-4,  are  naturally 
attenuated  under  certain  conditions.  Many  models  have  been  developed  to  simulate  the  fate 
and  transport  of  these  chemicals  in  different  media.  Unfortunately,  no  model  includes  every 
parameter  and  condition  and  is,  therefore,  unable  to  predict  exactly  what  will  happen  to  a 
certain  chemical  under  certain  conditions.  The  BIOPLUME  n™  model  was  developed  to 
predict  the  fate  and  transport  of  hydrocarbons  in  the  saturated  zone.  Results  from  this  study 
indicate  that  JP-4  is  naturally  attenuated  in  the  saturated  zone  and  certain  conditions  are 
favorable  to  the  process.  It  was  found  that  decreasing  hydraulic  conductivity,  increasing 
reaeration  rate  and  increasing  anaerobic  biodegradation  all  favor  natural  attenuation.  A 
general  mathematical  relationship  was  suggested  which  should  give  the  tield  engineer  a 
precursory  idea  whether  natural  attenuation  will  be  effective  in  reducing  contaminant 
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concentration  and  plume  travel  at  the  ate.  At  a  ate  wtere  these  mechanisms  can  be  verified, 
it  may  be  sufficient  to  rely  on  natural  attenuation  to  reaore  a  ate. 

Although  results  from  the  model  simulations  indicate  that  natural  attenuation  is 
enhanced  by  large  values  of  hydraulic  conductivi^,  reaeration,  and  anaerobic  biodegradation 
rates,  it  would  be  erroneous  to  state  that  under  these  circumstances,  natural  attenuation  should 
be  assumed  to  be  applicable  to  reduce  risk  at  a  site.  There  are  too  many  assumptions  built 
into  the  model,  too  many  simplifications  made,  and  too  many  parameters  excluded  to 
conclude  that  natural  attenuation  should  be  implemented  at  a  site  even  though  these  three 
parameters  indicate  a  favorable  attenuation  environment  All  the  factors  discussed  in  chapter 
two  must  be  considered  before  natural  attenuation  is  implemented. 

In  many  cases,  natural  attenuation  may  be  sufficient  as  a  remediation  option  and  where 
it  is  sufficient  it  ought  to  be  implemented.  Natural  attenuation  can  be  an  effective  and  cost 
efficient  method  for  site  restoration. 

Recommendations  for  Further  Research 

The  model  used  in  this  study  includes  numerous  parameters  that  were  not  investigated 
during  this  study.  Further  investigation  of  the  model  and  associated  parameters  may  prove 
useful.  The  code  for  the  model  is  available  and  it  may  be  interesting  to  improve  the  model 
by  adding  other  parameters  not  presently  included. 

Data  is  currently  being  gathered  at  three  Air  Force  sites  which  would  provide  real  data 
for  the  model.  An  investigation  of  the  capabilities  of  the  model  could  be  examined  using  this 
data. 
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In  light  of  the  current  budget  situation,  an  investigation  of  the  cost  savings  associated 
with  natural  attenuation  may  be  useful  for  justifying  use  of  this  alternative. 
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Appendix  A:  Glossary  of  Terms 

This  glossary  was  taken  from  Mobility  and  Degradarjf^n  nf  fh-yanic  Contaminants  in 
Subsurface  Environments  by  Warren  J.  Lyman  and  others. 

Anaerobe;  An  organian  that  does  not  require  air  or  free  oxygen  to 
maintain  its  life  process. 

Abiotic:  Referring  to  the  absence  of  living  organisms. 

Advection:  The  process  of  transfer  of  fluids  (vapors  or  liquids)  through  a  geologic  ftmnation 
in  response  to  a  pressure  gradient  which  may  be  caused  by  changes  in  barometric  pres^re, 
water  table  levels,  wind  fluctuations,  or  rainfall  percolation.  Advection  can  result  from  a 
thermal  gradient  caused  by  a  heat  source.  See  Darcy's  law. 

Aliphatic:  Of  or  pertaining  to  a  broad  category  of  carbon  compounds  distinguished  only  by  a 
straight,  or  branched,  open  chain  arrangement  of  the  consdment  carbon  atoms.  The 
carbon-carbon  bonds  may  be  saturated  or  unsaturated. 

Anisotropy:  The  dependence  of  property  upon  direction  of  measurement  (e.g.,  hydraulic 
conductivity,  porosity,  compressibility,  disperaon,  etc.). 

Aromatic:  Of  or  pertaining  to  organic  compounds  that  resemble  benzene  in  chemical 
behavior. 

Bentonite:  A  colloidal  clay,  largely  made  up  of  the  mineral  sodium 
montmorillonite,  a  hydrated  aluminum  silicate. 

Brownian  Movement:  Random  movement  of  molecules  or  colloids  su^nded  in  a  fluid. 

Biodegradation:  A  process  by  which  microbial  organisms  transform  or  alter  through 
en^matic  action  the  structure  of  chemicals  introduced  into  the  environment 

Biomass:  The  amount  of  living  matter  in  a  given  area  or  volume. 

Biota:  A  term  that  encompasses  the  spectrum  of  living  things  within  a  given  area. 

Biotic:  Of  or  pertaining  to  life  and  living  organisms.  Induced  by  the  actions  of  living 
organisms. 

Caliche:  A  layer  of  calcium  and  magnesium  carbonate  deposited  in  a 
near-surface  soil  horizon  by  evaporating  soil  water  or  groundwater.  It  may  occur  as  a  soft  thin 
soil  horizon,  as  a  hard  thick  bed  just  below  the  land  surface,  or  as  a  surface  layer  exposed 
by  erosion. 
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Capillarity:  The  pressure  difference  across  the  interface  of  two  immiscible  liquids  or  between 
a  liquid  and  a  solid.  In  a  porous  medium,  capillarity  is  zero  in  the  saturated  zone  and  less 
than  zero  in  the  vadose  zone.  See  surface  tenaon,  saturated  zone,  vadose  zone. 

Capillary  Fringe:  The  zoim  of  a  porous  medium  above  the  water  table  within  which  the 
porous  medium  is  saturated  but  is  at  less  than  atmospheric  pressure.  The  capillary  fringe  is 
considered  to  be  part  of  the  vadose  zone,  but  not  of  the  unsaturated  zone. 

Catalyst:  A  substance  that  alters  the  rate  of  a  ctemical  reaction  and  may  be  recovered, 
essentially  unaltered,  in  form  and  amount  at  the  end  of  the  reaction. 

Chemisorption:  A  process  in  which  weak  chemical  bonds  are  formed 
between  gas  or  liquid  molecules  and  a  solid  surface. 

Clay:  A  mineralogical  term  describing  a  family  of  aluminosilicate  minerals  formed  by  the 
decomporition  of  primary  minerals  (e.g.,  mica  and  feldspars)  and  re-composed  into  clay 
minerals.  A  textural  term  referring  to  particles  <  2  micrometers  in  diameter. 

Coalescence:  The  bonding  of  welded  materials  into  one  body,  or  the 
uniting  by  growth  in  one  body  as  particles,  gas  or  liquid. 

Colloid:  Particles  having  dimensions  of  10-10,000  angstroms  (1-1000  nanometers)  and  which 
are  dispersed  in  a  different  phase,  such  as  a  fluid  or  liquid. 

Complex:  A  combination  of  two  or  more  atoms  into  a  molecular  species,  usually  charged, 
and  existing  in  water  or  some  other  fluid. 

Compressibility:  The  change  in  volume  of  a  porous  medium  in  response  to  an  applied  stress 
which  is  counterbalanced  by  the  incompressibility  of  the  saturating  fluid  and  the  granular 
skeleton  of  the  porous  medium  which  it  saturams.  Compressibility  has  units  of  inverse  force 
(Pa-1). 

Constituents:  An  essential  part  or  component  of  a  ^stem  or  group: 
examples  are  an  ingredient  of  a  chemical  system,  or  a  component  of  an  alloy. 

Containment:  The  prevention  of  the  spreading  of  oil  or  other  hazardous  materials  by  the 
placing  of  booms  or  physical  barriers  and  the  use  of  absorbents,  gelling,  or  herding  agents  or 
other  materials  to  restrain,  entrap,  and  collect  a  spill. 

Corrective  Action:  The  removal  of  chemicals  and/or  contaminated 

soils,  objects  and  groundwater,  from  a  site,  or  other  clean-up  activities  designed  to  restore  the 
local  environment  to  acceptable  conditions.  Corrective  actions  may  include,  for  example, 
vacuum  extraction  of  the  vadose  zone,  soil  waging,  and  the  extraction  and  treatment  of 
contaminated  groundwater. 


Cytoplasm:  The  protoplasm  of  an  animal  or  plant  cell  external  to  the  nucleus. 


Darcy's  Law:  An  empirical  relation^ip  between  hydraulic  gradient  and  die  viscous  flow  of 
water  in  the  saturated  zone  of  a  porous  medium  untter  conditions  of  laminar  flow.  The  flux  of 
vapors  through  the  voids  of  the  vadose  zone  can  be  related  to  pressure  gradient  through 
the  air  permeability  by  Darcy's  law.  See  hydraulic  conductivity,  air  permeability,  hydraulic 
gradient,  pressure  gradient,  laminar  flow,  vadose  zone,  saturated  zone. 

Dehydrogenation:  Removal  of  hydrogen  from  a  compound. 

Density:  The  amount  of  mass  of  a  substance  per  unit  volume  of  the 
substance  (g/cm^). 

Desorption:  The  process  of  removing  a  sorbed  substance  by  the  reverse  of  adsorption  or 
absorption. 

Diffusion:  The  process  whereby  the  molecules  of  a  compound  in  a  single  phase  equilibrate  to 
a  zero  concentration  gradient  by  random  molecular  motion.  The  flux  of  molecules  is  from 
regions  of  high  concentration  to  low  concentration  and  is  governed  by  Pick’s  Second  Law. 

See  Pick's  Second  Law,  effective  diffusion  coefficient. 

Dispersion:  The  process  by  which  a  substance  or  chemical  spreads  and  dilutes  in  flowing 
groundwater  or  soil  gas.  On  a  microscale,  dispersion  is  due  to  mixing  within  individual  pores, 
mixing  between  pore  channels,  and  mixing  due  to  molecular  diffusion.  At  larger  scales, 
geologic  heterogeneity  and  anisotropy  cause  dispersion.  Dispersion  has  units  of  squared  length 
per  time  (cmVsec). 

Dissociation:  Separation  of  a  molecule  into  two  or  more  fragments  (atoms,  ions,  radicals)  by 
collision  with  a  second  body  or  by  the  absorption  of  electromagnetic  radiation. 

Dissolution:  Dissolving  of  a  material  in  a  liquid  solvent  (e.g.,  water). 

Dynamic  Viscosity:  The  measure  of  internal  friction  of  a  fluid  that  resists  shear  within  the 
fluid;  the  constant  of  proportionality  between  a  shear  stress  applied  to  liquid  and  the  rate  of 
angular  deformation  within  the  liquid,  having  units  of  mass  per  length  per  time  (gm/cm/sec). 

Effective  Di^usion  Coefficient:  The  constant  of  proportionality  in  Pick's  Second  Law  which 
is  dependent  on  tortuosity,  porosity,  and  moisture  content  and  properties  of  the  diffusing 
compound,  having  units  of  squared  length  per  time  (cmVsec).  See  tormosity,  porosity, 
moisture  content. 

Emulsification:  The  process  of  dispersing  one  liquid  in  a  second  immiscible  liquid. 
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Entrainment;  A  process  in  which  solid  particles  or  liquid  droplets  are,  by  force  of  friction 
with  a  passing  fluid  (e.g.,  air  or  water),  lifted  from  a  resting  place  and  carried  along  with  the 
flowing  fluid. 

Pick's  Second  Law:  An  equation  relating  the  change  of  concentration  with  time  due  to 
diffuaon  to  the  change  in  concentration  gradient  with  distance  from  the  source  of 
concentration.  See  diffusion,  effective  diffuaon  coefficient. 

Field  Capacity:  The  percentage  of  water  remaining  in  the  soil  2  or  3  days  after  gravity 
drainage  has  ceased  from  saturated  conditions. 

Fluid  Conductivity:  The  conaant  of  proportionality  in  Darcy's  law  relating  the  rate  of  flow 
of  a  fluid  through  a  cross-section  of  porous  medium  in  response  to  a  hydraulic  gradient.  Fluid 
conductivity  is  a  function  of  the  intrinsic  permeability  of  a  porous  medium  and  the 
kinematic  viscosity  of  the  fluid  which  flows  through  it.  Fluid  coiKlucdvity  has  units  of  length 
per  time  (cm/sec). 

Flux:  The  rate  of  movement  of  mass  through  a  unit  cross-sectional  area  per  unit  time  in 
response  to  a  concentration  gradient  or  some  advective  force,  having  units  of  mass  per  area 
per  time  (g/cm^*sec). 

Freundlich  Isotherm:  An  expression  relating  the  equilibrium  between  the  sorbed  chemical 
concentradon  and  its  aqueous  phase  concentration  through  a  power  law. 

Fugacity:  A  function  used  as  an  analog  of  the  partial  pressure  in  applying  thermodynamics  to 
real  ^sterns;  at  a  constant  temperature,  it  is  proportional  to  the  exponential  of  the  ratio  of 
the  chemical  potential  of  constituent  of  a  ^stem  divided  by  the  product  of  the  gas  constant 
and  the  temperature,  and  it  approaches  the  partial  pressure  as  the  total  pressure  of  the  gas 
approaches  zero. 

Fulvic  Acid:  A  term  of  varied  usage  but  usually  referring  to  the  mixture  of  organic 
substances  remaining  in  solution  upon  acidifrcation  of  a  dilute  alkali  extract  from  the  soil. 
Thus,  fulvic  acids  are  soluble  under  all  pH  conditions. 

Funicular  Zone:  A  narrow  band  above  the  groundwater  table  bounded 

above  by  the  capillary  rise  in  the  slimmest  continuous  pores,  and  bounded  below  by  the 

capillary  rise  in  the  widest  pores. 

Half-life;  The  time  required  for  half  of  a  substance  to  decay  or  alter  by  a  process,  (e.g., 
radioactive  decay,  biodegradation,  volatilization,  photolysis,  hydrolysis,  oxidation,  etc.). 

Henry's  Law:  The  relationship  between  the  partial  pressure  of  a  compound  and  the 
equilibrium  concentration  in  the  liquid  through  a  constant  of  proportionality  known  as  Henry's 
law  constant.  See  partial  pressure. 
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Heterogeneity;  The  dependence  of  [Hoperty  upcm  location  of  measurement  (e.g.,  hydraulic 
conductivity,  porosity,  compressibility,  di;q>ersion,  etc.)-  Heterogeneity  may  be  due  to  grain 
size  trends,  stratigraphic  contacts,  faults,  and  vertical  bedding. 

Homogeneity:  The  independence  of  property  widi  location  of  measurement  (e.g.,  hydraulic 
conductivity,  porosity,  compressibility,  (hq)ersion,  etc.). 

Humic  acid:  That  fraction  of  humic  substances  that  is  not  soluble  in  water  under  acid 
conditions  (below  pH  2),  but  becomes  solubte  at  greater  pH. 

Humic  substances:  A  general  category  of  naturally  occurring,  bit^enic  heterogeneous  mganic 
substances  that  can  generally  be  characterized  as  being  yellow  to  black  in  color,  of  high 
molecular  weight,  and  refractory.  Humic  substances  include  humin,  humic  acids  and  fulvic 
acids. 

Humin:  That  fraction  of  humic  substances  that  is  not  soluble  in  water  at  any  pH  value. 

Hydraulic  Conductivity:  The  constant  of  proportionality  in  Darcy's  law  relating  the  rate  of 
flow  of  water  through  a  cross-section  of  porous  medium  in  re^nse  to  a  hydraulic  gradient. 
Also  known  as  the  coefficient  of  permeability,  hydraulic  conductivity  is  a  function  of  the 
intrinsic  permeability  of  a  porous  medium  and  the  kinematic  viscosity  of  the  water  which 
flows  through  it.  Hydraulic  conductivity  has  units  of  length  per  time  (cm/sec). 

Hydraulic  Gradient:  The  change  in  piezometric  head  between  two  points  divided  by  the 
horizontal  distance  between  the  two  points,  having  dimensions  of  length  per  length  (cm/cm). 
See  piezometric  head. 

Hydraulic  Head:  A  measure  of  mechanical  energy  per  unit  weight  density  of  water  as  the 
sum  of  elevation  head  and  pressure  head,  having  units  of  length  (cm). 

Hydrocarbon:  One  of  a  very  large  group  of  ctemical  compounds  com¬ 
posed  only  of  carbon  and  hydrogen;  the  largest  source  of  hydrocarbons  is  from  petroleum 
crude  oil. 

Hydrolysis:  Decomposition  or  alteration  of  a  chemical  substance  by 
reaction  with  water. 

Hydrophobic:  Lacking  an  affinity  for,  repelling,  or  failure  to  absorb  or  dissolve  in,  water. 

Hysteresis:  The  dependence  of  the  state  of  a  system  on  direction  of  the  process  leading  to  it; 
a  non-unique  response  of  a  system  to  stress,  re^nding  differently  when  the  stress  is  released. 
Compressibility,  moisture  content,  soil  adsorption  and  unsaturated  hydraulic  conductivity 
exhibit  hysteretic  behavior. 
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niite:  A  group  of  2:1,  non-expanding,  hydrous  potassium-magnesium 
aluminosilicate  clay  minerals  of  intennediate  properties  between  kaolinite  and 
montmorillonite. 

Immersion:  Placement  into  or  within  a  fluid,  usually  water. 

Infiltration:  The  downward  movement  of  water  through  a  soil  from 

rainfall  or  from  the  application  of  artificial  recharge  in  response  to  gravity  and  capillaritj'. 

Intrinsic  Permeability:  A  measure  of  the  ease  with  which  a  porous  medium  transmits  a  fluid 
which  is  independent  of  the  properties  of  the  fluid  and  which  has  units  of  squared  length 
(cm2). 

Isotropy:  The  independence  of  a  property  with  direction  of  measurement  (e.g.,  hydraulic 
conductivity,  porosity,  compressibility,  dispersion,  etc.). 

Kaolinite:  A  1:1,  non-expanding,  hydrous  aluminosilicate  clay  mineral  formed  by  the 
alteration  v  ’  micas  and  feldspars. 

Kinematic  Viscosity:  The  ratio  of  the  dynamic  viscosity  of  a  fluid  to  its  density,  having  units 
of  squared  length  per  time  (cm2/sec). 

Locus  (pi,.  Loci):  In  this  report,  the  term  locus  is  used  to  refer  to  one  of  13  generic 
contaminated  environments  in  tlie  subsoil  region.  Each  locus  is  deHned  by  considering:  (I)  its 
position  relative  to  the  groundwater  table  (above,  on,  below);  (2)  the  phase  of  the  contaminant 
(liquid,  vapbr,  aqueous  solution,  sorbed  to  soil);  and  (3)  the  nature  of  the  local  natural 
matter  (unconsolidated  sediments,  fractured  rock,  flowing  water). 

Maeropore:  A  large  pore  in  a  porous  medium  which  may  be  formed  by 

physical  phenonema  or  biological  activity,  and  through  which  water,  or  other  fluids,  flows 

solely  under  the  influence  of  gravity,  unaffected  by  capillarity. 

Meniscus:  The  curved  surface  of  a  liquid  between  solid  boundaries  (e.g.,  capillary  tube, 
mineral  grains)  which  is  caused  by  the  surface  tension  of  the  liquid.  The  geometry  of  the 
meniscus  is  affected  by  the  surface  tension  of  the  liquid,  the  difference  in  density  between 
the  liquid  and  overlying  air,  the  distance  between  the  solid  boundaries,  and  the  hydrophobic 
nature  of  the  liquid. 

Mica:  A  family  of  platy,  igneous  and  metamorphic,  aluminosilicate  minerals  which  weather 
and  form  clays.  Micas  separate  readily  into  thin  sheets  or  flakes. 

Microbe:  A  microorganism,  especially  a  bacterium. 
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Microorganisms:  Microscopic  organisms  including  bacteria,  protoioans,  yeast,  fungi,  viruses 
and  algae. 

Mobilization;  The  process  or  processes  by  which  a  liquid  contaminant  in  a  locus  is  made 
more  mobile  in  the  locus  or  more  transferable  between  loci. 

Moisture  Content:  The  amount  of  water  lost  from  the  soil  upon  drying  to  a  constant  weight, 
expressed  as  the  weight  per  unit  weight  of  dry  soil  or  as  the  volume  of  water  per  unit  bulk 
volume  of  the  soil.  For  a  fully  saturated  medium,  moisture  content  equals  the  porosity;  in 
the  v  adose  zone,  moisture  content  ranges  between  zero  and  the  porosity  value  for  the  medium. 
See  porosity,  vadose  zone,  saturated  zone. 

Mole  Fraction:  The  ratio  of  the  number  of  moles  of  a  substance  in  a  mixture  or  solution  to 
the  total  number  of  moles  of  all  the  components  in  the  mixture  or  solution. 

Monovalent:  A  radical  or  atom  whose  valency  is  1. 

Montmorillonite:  A  2:1,  expanding,  hydrous  magnesium  aluminosilicate  clay  mineral 
exhibiting  pronounced  swelling-r''rinkage  behavior  and  high  plasticity  and  cohesion. 

Oxidation:  A  chemical  reaction  thiiv  increases  the  oxygen  content  of  a  compound,  or  raises 
the  oxidation  state  of  an  element. 

Oxidation  Potential;  The  difference  in  potential  between  an  atom  or  ion  and  the  state  in 
which  an  electron  has  been  removed  to  an  infinite  distance  from  this  atom  or  ion. 

Partial  Pressure:  The  portion  of  total  vapor  pressure  in  a  ^stem  due  to  one  or  more 
constituents  in  the  vapor  mixture. 

Percolation,  Soil  Water:  The  downward  movement  of  water  through 

soil.  Especially,  the  downward  flow  of  water  in  saturated  or  nearly  saturated  soil  at  hydraulic 

gradients  of  the  order  of  1.0  or  less. 

Phreatic  Level:  The  groundwater  level  that  would  be  seen  in  an  observation  well  projecting 
down  into  an  aquifer. 

Polymerization;  The  bonding  of  two  or  more  monomers  to  produce  a 
polymer. 

Porosity:  The  volume  fraction  of  a  rock  or  unconsolidated  sediment  not  occupied  by  solid 
material  but  usually  occupied  by  water  and/or  air.  Porosity  is  a  dimensionless  quantity. 


73 


Redox:  A  chemical  reaction  in  which  an  atom  or  molecule  lossea/gains  electrons  uVfrom 
another  atom  or  molecule.  Also  called  oxidation-reduction.  Oxidation  is  the  loss  of  electrons; 
reduction  is  tte  gain  in  electrons. 

Reduction:  A  chemical  reaction  in  which  an  atom  of  molecule  gains 
electrons.  Sometimes  results  by  reaction  of  the  substance  with 
hydrogen. 


Refractory:  A  nonspecific  characteristic  of  some  chemicals  implying  resistance  to 
biodegradation  or  other  degradation  or  treatment  processes. 

Residence  Time:  The  average  time  that  water  remains  in  a  porous  medium  or  a  particle 
remains  in  a  reservoir.  Residence  time  is  calculated  as  the  ratio  of  reservoir  volume  to  total 
water  inflow  rate  having  units  of  time  (sec). 

Residual  Saturation:  The  amount  of  water  or  oil  remaining  in  the  voids  of  a  porous 
medium  and  held  in  an  immobile  state  by  capillarity  and  dead-end  pores. 

Salinity:  The  quantity  of  anions  and  cations  in  water,  usually  between  33  and  37  parts  per 
thousand  in  sea  water. 

Sinter:  A  chemical  sedimentary  rock  deposited  by  precipitation  from 
mineral  waters,  especially  siliceous  sinter  and  calcareous  sinter. 

Slurry:  A  thick  mixture  of  liquid,  especially  water,  and  any  of  several 
frnely  divided  substances,  such  as  cement  or  clay  particles. 

Solubility;  The  maximum  amount  of  mass  of  a  compound  that  will 

dissolve  into  a  unit  volume  of  solvent,  usually  water,  having  units  of  mass  per  volume 

(gm/cm3). 

Sorption:  A  general  term  used  to  encompass  the  process  of  absorption, 
adsorption,  ion  exchange,  and  chemisorption. 

Specific  Gravity:  The  ratio  of  the  weight  of  a  given  volume  of  the 
material  at  4°C  (or  some  stated  temperature)  to  the  weight  of  an  equal  volume  of  distilled 
water.  Materials  with  qiecifrc  gravity  of  less  than  1  will  float  on  water,  materials  with 
specific  gravity  over  1  will  sink  in  water. 

Surface  Tension:  A  measure  of  the  interfacial  tension  due  to  molecular 

attraction  between  two  fluids  in  contact  ot  between  a  liquid  in  contact  with  a  solid,  having 

units  of  mass  per  squared  time  (dyne/cm). 
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Surfactant:  Also  Surface  Active  Agent  A  chemical  material  which  pro¬ 
vides  a  "linking  action"  between  two  materials,  such  as  oil  and  water,  which  normally  resist 
mixing  or  readily  joining  in  solution.  Surfactants  provide  die  emulaficadon  forces  which 
allow  oil  and  water  to  mix  and  remain  in  either  oil-in-water  or  water-in-oil  solutions. 

Tortuosity;  The  ratio  of  path  length  through  a  porous  medium  to  the 
straight-line  flow  path  which  describes  the  geometry  of  the  porous 
medium.  Tortuosity  is  a  dimensionless  parameter  which  ranges  in 
value  from  1  to  2. 

Transpiration:  The  release  of  water  withdrawn  from  the  soil  by  plants 
during  photo^nthesis  and  other  life  processes. 

Unsaturated  Zone:  The  portion  of  a  porous  medium,  usually  above  the 
water  table  in  an  unconfined  aquifer,  within  which  the  moisture  con¬ 
tent  is  less  than  saturation  and  the  capillary  pressure  is  less  than 
atmospheric  pressure.  The  unsaturated  zone  does  not  include  the 
capillary  fringe. 

Vadose  Zone:  The  portion  of  a  porous  medium  above  the  water  table 
within  which  the  capillary  pressure  is  less  than  atmospheric  and  the 
moisture  content  is  usually  less  than  saturation.  The  vadose  zone 
includes  the  capillary  fringe. 

Vapor  Pressure;  The  equilibrium  pressure  exerted  on  the  atmosphere  by 
a  liquid  or  solid  at  a  given  temperature.  Also  a  measure  of  a  sub¬ 
stance's  propensity  to  evaporate  or  give  ofr  flammable  vapors.  The 
higher  the  vapor  pressure,  the  more  volatile  the  substance. 

Vermiculite:  A  2:1,  hydrous  ferro-magnesium  aluminosilicate  clay  mineral  similar  in  structure 
and  property  to  montmorillonite. 

Vitrification:  Formation  of  a  glassy  or  noncrystalline  material. 

Volatilization:  The  process  of  transfer  of  a  chemical  from  the  water  or 
liquid  phase  to  the  air  phase.  Solubility,  molecular  weight,  and  vapor  pressure  of  the  liquid 
and  the  nature  of  the  air-liquid/water  interface  affect  the  rate  of  volatilization.  See  solubility, 
vapor  pressure. 

Water  Table:  The  water  surface  in  an  unconfined  aquifer  at  which  the 
fluid  pressure  in  the  voids  is  at  atmospheric  pressure. 

Wilting  Point:  The  point  at  which  a  plant  wilts,  no  longer  able  to  withdraw  water  from  a  soil 
or  sediment  to  support  transpiration  processes  and  retain  turgor  pressure. 
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Appendix  B 
Sample  Input  Data 
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Bmmplu  Input  Data  -  IUX1IM1.  MOOPaMMltlOW 


Caxd  1 


1  1 

0  0 

11 

0 

20  2000 

1 

1 

7 

1  100 

1  9  X  0 

0 

c.oo 

.001 

0.3  10 

0 

0 

0  50 

50  0.1  0.5 

1.0 

0  0  0  0  0 

5  10 

6  6  -.0002  100  0 

0  .0025 

0  25.0 


-  Cacd  2 

-  Cacd  3 

-  Caxd  4 

-  Data  dat  1 

-  Data  dat  2 

-  Data  dat  3 

-  Data  dat  4 


0  0.0 

1  1.0 

00000000000 
01111111110 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
01111111110 
00000000000 


1 

1 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


1.0  0.0  8.0 
1.0 

0.  0.  0.  0. 

100.  100.  100.  100. 


97.  97.  97.  97. 

0.  0.  0.  0. 


0.0 


0. 

100. 


97. 

0. 


0 


0. 

100. 


97. 

0. 


0.  0.  0. 

100.  100.  100. 


97.  97.  97. 

0.  0.  0. 


0. 

0. 


0. 

0. 


-  Data  dat  5 

-  dtart  of  Data  dat  C 


-  Bnd  of  Data  dat  ( 

-  Data  dat  7 

-  Start  of  Data  dat  8 


-  dnd  of  Data  dat  8 


0  0.0 

0  8.0 


-  Data  Sat  9 

-  Data  dat  10 
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Appendix  C 
Sample  OuQKit  Data 


Appendix  C:  Sample  Output  Data 
1  BIOPLUME  n 

ICONTAMINANT  TRANSPORT  UNDER  THE  INFLUENCE  OF  OXYGEN  LIMITED  BIODEGRADATION 

0  Bioplume  II  -  Sample  Data 

0  INPUT  DATA 

0  GRID  DESCRIPTORS 

(NUMBER  OF  COLUMNS)  =  11 

(NUMBER  OF  ROWS)  =  20 

EL  (X-DISTANCE  IN  FEET)  =  50.0 

YDEL  (Y-DISTANCE  IN  FEET)  =  5OJ0 
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TIME  PARAMBICRS 


NTIM  (MAX.  NO.  (XF  TIME  STEPS)  »  1 

NPMP  (NO.  C»>  PUMPING  PERIODS)  »  1 

PINT  (PUMPING  PERIOD  IN  YEARS)  6j000 

TIMX  (TIME  INCREMENT  MULTIPLIER)  »  EO 

TINTT  (INITIAL  TIME  STEP  IN  SEC.)  =  0. 

0  HYDROUXJIC  AND  CHEMICAL  PARAMETERS 

S  (STORAGE  COEFFICIENT) 

POROS  (EFFECTIVE  POROSITY)  =  300 


BETA  (LONGITUDINAL  DISraRSIVITY)  =  IOjO 
DLTRAT  (RATIO  OF  TRANSVERSE  TO 
LONGITUDINAL  DISPERSIVTTY) 

ANFCTR  (RATIO  OF  T-  YY  TO  T-XX) 

0  EXECUTION  PARAMETERS 


NITP  (NO.  OF  ITERATTON  PARAMETERS)  =  7 

TOL  (CONVERGENCE  CRITERIA  -  ADD*)  =  0010 

ITMAX  (MAX.  NO.  OF  ITERATIONS  -  ADD*)  =  100 
CELDIS  (MAX.  CELL  DISTANCE  PER  MOVE 
OF  PARTICLES  •  M.O.C.) 

NPMAX  (MAX.  NO.  OF  PARTICLES)  =  2000 

NPTPND  (NO.  PARTICLES  PER  NODE)  =  9 


1 

0  PROGRAM  OPTIONS 


NPNT  (TIME  STEP  INTERVAL  FOR 
COMPLETE  PRINTOUT) 

NPNTMV  (MOVE  INTERVAL  FOR  CHEM. 

CONCENTRATION  PRINTOUT) 

NPNTVL  (PRINT  OPTION-VELOCITY 
0=NO;  1=FIRST  TIME  STEP, 

2=ALL  TIME  STEPS)  =  0 

NPNTD  (PRINT  OPnON-DISP.COEP. 

0=NO;  1=FIRST  TIME  STEP 
2=ALL  TIME  STEPS) 

NUMOBS  (NO.  OF  OBSERVATION  WELLS 
FOR  HYDROGRAPH  PRINTOUT)  *  1 


jOOOOOO 

.10 

IBOOOOO 


300 


1 

0 

0 
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NREC  (NO.  OF  PUMPING  WELLS)  =  1 

NCOOES  (FOR  NOLffi  IDENT.)  >  1 

NPNCHV  (PUNCH  VEL(X:iTY)  =  0 

NPDELC  (PRINT  OPT.-CONC.  CHANGE)  =  0 

0  REACTION  TERMS 


DK  (DISTRIBUTION  COEPFICIEhrn 
RHOB  (BULK  DENSrrV  OF  SOLIDS) 

RF  (RETARDATION  FACTOR) 

THALF  (HALF  LIFE  OF  raCAYJN  SEC)  = 
DECAY  (DECAY  CONSTANT=LN  2/THALF) 


s  jOOO(X£400 
=  joooo(£+oo 

iOOOOEtOl 

jOOOOCffi+OO 

s  jOOOOOEfOO 


0 


DECAY  TERMS 


DECl  (ANAEROBIC  DECAY  COEFF.  )  =  j0000(»+OO 

DEC2  (REAERATTON  DECAY  COEFF.)  =  J00000E4O0 


1  STEADY-STATE  FLOW 

TIME  INTERVAL  (IN  SEC)  FOR  SOLUTE-TRANSPORT  SIMULATION  *  .18935E+09 

0  LOCATION  OF  OBSERVATION  WELLS 

NO.  X  Y 
1  6  10 

0  LOCATION  OF  PUMPING  WELLS 

X  Y  RATEdN  CFS)  C(»IC.  CONC(02) 

6  6  -20E-03  lOOBO  .00 

0  AREA  OF  ONE  CELL  =  2500. 
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0 


X-Y  SPACINO: 
SOjOOO 
SOXXX) 


ITRANSMISSIVITY  MAP  (FT*FT/SEC) 

OXXH&iOO  OjOOEfOO  OJXB+00  0.0(»^  0j0(»+O0  OjO(»+00  OJ)OBHX)  OXXB+OO  OXMBfOO 


0X)(ffi-KI0 

Z5<m4i3 

2.50E-O3 

250E-03 

2J0E-03 

2Jm^ 

2J0B-03 

2J()E-03 

OJOOB*00 

0X)QE+00 

2J(m-03 

2J(MB-03 

2J0E-03 

2.5(£-03 

2J(m-C3 

2.50E^ 

2.soe-o3 

0J)08+00 

0.0(ffi400 

230B-03 

2.50E-03 

2.50E-03 

2.5(m^ 

2308-03 

230EM 

2.5(»^ 

0008*00 

0.0QEf00 

2.50E^ 

2J0E-03 

230E^ 

2J0E-O3 

23(^03 

230E4a 

23084)3 

0008*00 

0.0(HB+00 

23(ffi-03 

2J0E-03 

2JCEr03 

2J0E-O3 

230E-03 

23084)3 

2.5flE-03 

0008+00 

0X)QE+00 

2J0E-03 

2j(ffi-03 

2.50E-03 

23(HB-03 

230Er€3 

23084)3 

2.50e^ 

0008*00 

0X)0E+00 

2.50E-03 

23(ffi-03 

2JQE-03 

2J(£<03 

230E433 

23CB4)3 

230e-03 

0008*00 

O.OOE+00 

2.50E-03 

2J(ffi-03 

2.50E-03 

2JHffi-03 

23(SrQ3 

23084)3 

230E^ 

0008*00 

OjOOE+00 

230E-03 

2J(m-03 

2J(m-03 

2J<m^ 

230E4>3 

2J(£-03 

230E-03 

0008*00 

O.OOE+00 

2.50E-03 

2J0E^ 

2J<m-03 

2ii(«-03 

230Em 

2J<B-03 

2J(»-03 

0008+00 

O.OOE-IOO 

2.50E-03 

2J0E-03 

2.50E-03 

2J0E^ 

230E-03 

25(£-03 

2.5(SB-03 

0008+00 

O.OOE+00 

230E-03 

2J0E-03 

2J0E-03 

230E^ 

230E-03 

250E^ 

2J0E-03 

oo(m*oo 

o.o(m+oo 

2J0E-03 

2J(m-03 

2.50E^ 

2.50E-03 

2.5(»-03 

2J0E-03 

23(£-03 

0008+00 

O.OOE+00 

2.50E-03 

230E-03 

2J(ffi-03 

2.50E-03 

2.5(B-03 

23084)3 

23(m-03 

0.008*00 

OX)OE+O0 

2J0E-03 

2.5(m-03 

2J0E-03 

230E-03 

2.5(ffi-03 

23iX4>3 

2.50E-03 

0008+00 

0.0(»-tXX) 

2.50E-03 

2.5QE-03 

2.50E-03 

2.50E-03 

25(B-03 

250E-03 

2306-03 

o.o(m+oo 

O.OOE+00 

250E-03 

2J0E-03 

2J0E-03 

2.50E^3 

2.50E-03 

2J0E-03 

23(JE-03 

0.008+00 

O.OOE+00 

2J0E-03 

2J0E-O3 

23(ffi-03 

2.5(»-03 

2.50E-03 

250E^G 

23084)3 

o.0(m*oo 

0.0OE-K)0 

O.OOE+00 

0X)0Ef00 

O.OOB+00  OJOOE+00  0X)0Ef00  O.OOE400  OX)OE+O0  OXXB+OO 

lAQUIFER  THICKNESS  (FT) 


.0  .0  .0  .0  .0  X)  .0  jO  D 

.0  25X)  25.0  25X)  25.0  25j0  25j0  25j0  25.0  25X)  X) 

X)  25.0  25.0  25.0  25.0  25X)  25.0  25.0  25.0  25.0  X) 

.0  25.0  25.0  25.0  25X)  25.0  25.0  25.0  25.0  25.0  X) 

.0  25X)  25.0  25.0  25.0  25X)  25X)  25.0  25.0  25.0  X) 

.0  25.0  25.0  25.0  25.0  25.0  25X)  25.0  25.0  25.0  X) 

.0  25.0  25.0  25.0  25.0  25.0  25X)  25.0  25X)  25X)  .0 

X)  25.0  25.0  25.0  25X)  25.0  25.0  25X)  25.0  25X)  X) 

.0  25.0  25X)  25.0  25.0  25X)  25X)  25X)  25.0  25.0  .0 

.0  25.0  25.0  25X)  25.0  25.0  25.0  25.0  25X)  25X)  X) 

.0  25.0  25.0  25.0  25.0  25.0  25.0  25X)  25.0  25X)  X) 

.0  25X)  25.0  25X)  25.0  25X)  25X)  25X)  25X)  25X)  X) 

.0  25X)  25.0  25X)  25X)  25X)  25.0  25X)  25X)  25.0  X) 

.0  25.0  25X)  25X)  25X)  25.0  25.0  25.0  25X)  25.0  .0 


X)  X) 
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jO  2SjO  2SjO  2SJ0  2SjO  2SjO  2SjO  2SjO  25j0  2SX)  jO 

X)  2SjO  2SjO  25j0  25j0  25j0  25j0  25j0  25j0  2SjO  jO 

jO  25j0  25j0  25j0  25i>  25j0  25j0  25j0  25j0  25X)  jO 

X)  25X)  25X)  25i)  25j0  25j0  25J)  25i)  25X)  25X)  X) 

X)  25X)  25X)  25X)  25X)  25X)  25X)  25X)  25X)  25X)  X) 

X).O.OX)X)X)X)X)X)X)X) 


IDIFPUSE  RECHARGE  AND  DISCHARGE  (PT/SEC) 

OXXBfOO  OXXB+OO  OXKB+00  OXWEfOO  0XK£+00  OXX^-HW  OjOOB+00  OjOOB*€0  ODOE+OO  OXKB+OO 
OJOCB*CO 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

O.OOE-HX)  0.0(£400  OXWE^  OXXm+OO  OXX)E-K»  O.OOE400  0J00iE*O0  OXKB+OO  O.OOB+00  O.OOE-MX) 
O-OOE-HX) 


IPERMEABILTY  MAP  (FT/SEC) 

O.OOE-HX)  O.OOE+00  OXME+OO  0.00E-K)0  OXXB+OO  OX)OEtOO  0.0(ffi400  OXX^-HX)  0XX)E+00  0X)0E+00 
0.0<^+00 


O.OOE+00 

l.OOE-04 

l.OOE-04 

lX)0E-04 

l.OOE-04 

1.0(«E-04 

lX)0E-04 

l.OOE-04 

l.OOE-04 

lX)0E-04 

OXKE+00 

O.OOE+00 

OXX)E-fOO 

lX)0E-04 

l.OOE-04 

1X)OE-04 

l.OOE-04 

l.OOE-04 

1XKB04 

lX)0E-04 

lX)0E-04 

1.0(m-04 

O.OOEfOO 

OXKB+OO 

l.OaE-04 

l.OOE-04 

l.OOE-04 

lX)0E-04 

l.OOE-04 

lX)0E-04 

1.0(£-04 

1.0(£-04 

lX)(ffi-04 

0.0(£^ 

0XX)E+00 

l.OOE-04 

l.OOE-04 

l.OOE-04 

lX)0E-04 

l.OOE-04 

1.00E^ 

lX)0E-04 

\XKB4yi 

lX)(B-04 

OXXBfOO 

lX)0E-04 

l.OOE-04 

l.(X)E-04 

1X)OE-04 

lX)0E-04 

1X)OE-04 

lX)0E-04 

lX)0E-04 

lXXH-04 

OXKB+OO 
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OXME-KX) 

O.0OE-i4)O 

1.0(B-04 

lJOOE-04 

IjOOE-04 

lBQE-04 

lBQB-04 

lBQE-04 

lBQE-04 

1jO(£-04 

ljm£-04 

0.0QE-f00 

OBOE-tOO 

lBOE-04 

1.0(£-04 

lj0(ffi-04 

lBOE-04 

1.0(»-04 

lB(£-04 

1jO(£-04 

1jO(£-04 

lB(£-04 

O.OOE-fOO 

O.OOE+00 

1.0(£-04 

1.0(^-04 

1jOOE-04 

l.OOE-04 

lBQE-04 

lB(£-04 

1jO(£-04 

lBQE-04 

lB(£-04 

0.00E-K)0 

O.OOE400 

l.OOE-04 

l.OOE-04 

1.0(»E-04 

IjOOE-04 

1.0(ffi-04 

lB(B-04 

1.0(m-04 

1jO(£-04 

1jO(»B-04 

O.OOE-iBO 

O.OOEfOO 

lBOE-04 

l.OOE-04 

1B0E-04 

IjOOE-04 

lB(»-04 

l.OOE-04 

lB(£-04 

1jO(£-04 

1jO(£-04 

O.OOE+00 

O.OOE+00 

lBOE-04 

1.00E-04 

l.OOE-04 

lBOE-04 

lBOE-04 

1.0(£-04 

lB(»E-04 

1jO(B-04 

1jO(B-04 

O.OQE-fOO 

O.OOE+00 

l.OOE-04 

l.OOE-04 

l.OOE-04 

l.OOE-04 

lJOOE-04 

lBOE-04 

1jO(KE-04 

1.0(£-04 

lB(K-04 

0.00E-K)0 

O.OOE+00 

l.OOE-04 

l.OOE-04 

lBOE-04 

l.OOE-04 

lBOE-04 

l.OOE-04 

l.OOE-04 

1.0(®-04 

1.0(£-04 

0.00E-K)0 

O.OOE+00 

l.OOE-04 

1.00E-04 

1.00E-04 

l.OOE-04 

lBOE-04 

1.00E-04 

1.0(ffi-04 

1.0(B-04 

IjOOE-04 

0.00E4B0 

0.00E-»<X) 

l.OOE-04 

l.OOE-04 

l.OOE-04 

1X)OE-04 

l.OOE-04 

l.OOE-04 

1.0(£-04 

lB(£-04 

ixxm-04 

O.OOE+00 

O.OOE+00 

l.OOE-04 

l.OOE-04 

l.OOE-04 

l.OOE-04 

lBOE-04 

lBOE-04 

lBOE-04 

l.OOE-04 

lBOE-04 

O.OOE+00 

O.OOE+00 

l.OOE-04 

1.00E-04 

l.OOE-04 

l.OOE-04 

l.OOE-04 

lBOE-04 

lBOE-04 

l.OOE-04 

1B0E-04 

O.OOE+00 

l.OOE-04 

1.0(ffi-04 

1.00E-04 

1.00E-04 

1.00E-04 

lB(m-04 

lBOE-04 

1jO(£-04 

lB(B-04 

O.OOE+00 

O.OOE+00  O.OOE+00  O.OOE+00  O.OOE+00  O.OOE+00  O.OOE+00  O.OOE+00  O.OOE+00  O-OOE+OO  O.OOE+00 
O.OOE+00 

0  NO.  OF  FINITE-DIFFERENCE  CELLS  IN  AQUIFER  =  162 

AREA  OF  AQUIFER  IN  MODEL  =  .4050(®+{)6  SQ.  FT. 

NZCRIT  (MAX.  NO.  OF  CELLS  THAT  CAN  BE  VOID  OF 

PARTICLES;  IF  EXCEEDED.  PARTICLES  ARE  REGENERATED)  =  3 


INODE  IDENTIFICATION  MAP 

00000000000 
0  1111111110 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
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00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
0  1111111110 
00000000000 


0  NO.  OF  NODE  IDENT.  CODES  SPECIFIED  «  1 
0  THE  FOLLOWING  ASSIGNMENTS  HAVE  BEEN  MADE: 

CODE  NO.  LEAKANCE  SOURCE  CONC.  02CONC  RECHARGE 
0  1  .lOQE+01  XX)  8jOO 


IVERTICAL  PERMEABIUTY/THICKNESS  (FT/(FT*SEC)) 


0.00Ef00 

0.0OE-M)0 

O.OOE+00 

O.OOE+00 

O.OOE+00 

O.OOE-HOO 

O.OOB+00 

0.00E+00 

O.OOE+00 

O.OOE+00 

0.00E+00 

O.OOE+00 


0.00E^ 

l.OOE+00 

O.OOE+00 

O.OOE+00 

0.00E+00 

O.OOE+00 


0X)0E+00 

l.OOE+00 

0.00E400 

0.00E-f00 

OJOOE+00 

O.OOE+00 


O.OOE+00 

l.OOE+00 

oxm+oo 

ODOE-^00 

OJOGE+OO 

O.OOE+00 


O.OOE+00 

1.0(®+00 

OjOOE+OO 

0.0(ffi+00 

O.OOE+00 

O.OOE+00 


0.0(£+00  OJXS-KX) 
IDOE+OO  IDOE+OO 

o.ooE+00  o.ooe+oo 

0.00E-HX)  O.OOE+00 
O.OOE-fOO  0J)0E+00 
O.OOE+00  0.00E+00 


0D(S+00 

1j0(B+00 

ODOE+00 

0.0(£-H)0 

O.OQE+00 

O.OOE+00 


0.0(ME+00 

1.0(ME-f00 

OD(£400 

OBOEfOO 

0XX£-K)0 

0.0(«B400 


O.OOE+00 

IjOOE-HX) 

OBOE+OO 

0.0(®+0C 

OXX>E400 

O.OOE+00 
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0iX)E-KX) 

0i)0E+00 

p 

1 

0j0QE400 

OjOOE+00 

ojooe-foo 

ojxm^ 

OjO(£400 

OJOQE400 

OJ0QE4OO 

1 

o 

0.00E-K)0 

O.OOE-fOO 

1 

o 

0.00E-f00 

oxxm-hoo 

OjOGE-KX) 

osm+oo 

ojom+oo 

0X)0E400 

0.0(£400 

OJOQE400 

0.0(m+00 

0.0(^-K)0 

O.OOE+00 

1 

o 

0.0(£-f00 

O.OOE+00 

0.0(£400 

0iX»+00 

OjOOE400 

OB(£400 

0.00E400 

0.0OE-K)0 

O.OOE-KOO 

Oj0(ffi-K)0 

OjO(»+00 

OXXffi+OO 

ODOE-KX) 

OiXS+OO 

OX)QE400 

OX)OE400 

0XXS400 

OjO(£400 

O.OOE-KOO 

0.00E4O0 

0.0QE-K)0 

0.0(^-KX) 

OjOOE+OO 

O.OOE+00 

OBCB+OO 

0D(£400 

0X)(£400 

OjOOE400 

0.0(£400 

0.0OE4O0 

0.00E-K)0 

o 

1 

1 

o 

0.0(£-f00 

0X)0E+00 

0.0(£-f00 

O.0OE-K)O 

0.0(^400 

0j0Oc4O0 

O.0QE4O0 

O.OOE+00 

O.OOE+00 

O.OOE+00 

0.00E400 

1 

o 

1 

1 

o 

OOOE+OO 

0.0(^400 

0.0(£400 

0.0(B400 

0.00E4O0 

0.0OE4O0 

O.OOE+00 

O.OOE+00 

0.0(£-»00 

OjOOE+OO 

0.00E-KX) 

0.0(^400 

I 

1 

0.0(£400 

0.0OE4O0 

0X)(£400 

O.OOE+00 

0.00E4O0 

O.OOE+00 

O.0OE4O0 

0J0€E■^O0 

O.OOE+00 

I 

1 

OJ)OE400 

0.0(B400 

0.0OE4O0 

o 

1 

O.OOE+00 

O.OOE+00 

1 

o 

1 

o 

O.OOE+00 

O.OOE+00 

1 

o 

0.00E400 

0.00E400 

o 

1 

O.0OE4O0 

O.OOE+00 

O.OOE+00 

0.00E4<X) 

1 

o 

1 

o 

O.OOE+00 

O.OOE+00 

0.00E400 

O.0OE4O0 

0.00E400 

0.00E400 

O.OOE+00 

0.00E-f00 

O.OOE+00 

O.OOE+00 

O.OOE+00 

O.OOE-^00 

o 

I 

0.0OE40O 

o 

1 

I 

o 

0.00E400 

O.OOE+00 

O.OOE+00 

l.OOE+00 

l.OOE+00 

l.OOE+00 

l.OOE+00 

l.OOE+00 

1.00E400 

1.00E400 

1X)0E400 

1.00E400 

O.OOE+00 

O.OOE+00 

1 

o 

O.0OE4OO 

1 

o 

O.OOE+00 

O.OOE+00 

0.00E400 

0.00E400 

0.00E400 

0.00E400 

IWATER  TABLE 

0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0. 

0.  100.  100.  100.  100. 100.  100.  100.  100.  100.  0. 

0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0. 

0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0. 

0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0. 

0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0. 
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il.OOEUM) 
O.OOIi  i  (K) 

O.dOE !  (H) 

O.OOEKK) 

O.OOE+{K) 

O.OOE+OO 

O.OOEtOO 

p 

o 

rr 

-r 

8 

O.OOE+OO 

O.OOE+00 

0.00E+00 

O.dOF, !  (X) 
0.0()!v((K) 

0.()0E.i(K) 

O.OOf<+(X) 

O.OOE.+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

0.00E+00 

0.00E+00 

O.OOE+00 

().()0!vf(K) 

o.oon^oo 

().()OE+(K) 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+00 

O.OOE+00 

O.OOE+00 

().0()E+()0 

O.OOhfOO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

0.00E+00 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+00 

O.OOE+00 

O.OOEfOO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+Ob 

0.00E+00 

O.OOE+00 

0.00E+00 

O.OOE+00 

O.OOE+00 

O.OOE+00 

O.OOE+00 

O.OOE+OO 

O.OOE+00 

0.00E+00 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+00 

O.OOE+00 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

0.00E+00 

O.OOE+OO 

0.00E+00 

O.OOE+OO 

0.00E+00 

O.OOE+OO 

O.OOE+00 

O.OOE+00 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

0.00E+00 

O.OOE+00 

O.OOE+OO 

O.OOE+00 

O.OOE+00 

O.OOE+00 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+0O 

O.OOE+OO 

O.OOE+00 

O.OOE+00 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+00 

O.OOE+OO 

O.OOE+00 

O.OOE400 

O.OOE+00 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+00 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+00 

O.OOE+00 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

0.00E+00 

O.OOE+OO 

O.OOE+OO 

O.OOE+00 

0.00E+00 

0.00E+00 

O.OOE+00 

O.OOE+00 

O.OOE+OO 

O.OOE+OO 

l.OOE+00 

1.00E+00 

l.OOE+OO 

l.OOE+00 

1.00E+00 

1.00E+00 

l.OOE+00 

1.00E+00 

l.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

0.00E+00 

O.OOE+OO 

O.OOE+00 

0.00E+00 

O.OOE+00 

1  WATER  TABLE 


0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

100. 

100 

.  100 

1.  100. 100.  100. 

100, 

,  100.  100 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0.. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 
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i). 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

{). 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

C, 

0. 

0. 

0. 

0. 

0. 

« 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

’0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

97. 

97. 

97. 

97. 

97.  97.  97.  ' 

97. 

97. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

IHEAD  DISTRIBUTION  -  ROW 
NUMBER  OF  TIME  STEPS  =  0 

TIME(SECONDS)  =  .00000 
TIME(DAYS)  =  .OOOOOE+00 
'nME(YEARS)  =  .OOOOOE+00 


0  .0000000  .0000000  .0000000  .0000000  .0000000  .0000000 
.00000000  0 


1 1:* 


0  100.0000000  100.0000000  100.0000000  100.0000000  100.0000000  100.0000000  100. 
100.0000000  100.00000000  0 
0 
0 
0 


.0000000 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0  97.0000000  97.0000000  97.0000000  97.0000000  97.0000000  97.0000000  97.0000000  97X)000000 
97.0000000  0 

0  .0000000  .0000000  .0000000  XXXXXXW  JOOOOOOO  .0000000  .0000000  DOOOOOO  .0000000 

0 


IFTERATION  PARAMETERS 
.6168S0E-02 
.14404(^-01 
J36346E-01 
.785398E-01 
.193397 
.428249 
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1.00000 


.000000 

.000000 

.000000 

.000000 

.000000 

.000000 

.000000 

.000000 

.000000 

.000000 

.000000 

.000000 

.000000 


ICONCENTRATION  OF  CONTAMINANT 
NUMBER  OF  TIME  STEPS  =  0 

TlME(SECONDS)  =  .00000 

CHEM.TIME(SECONDS)  =  .OOOOOE+00 
CHEM.TIME(DAYS)  =  .OOOOOE+00 
TIME(YEARS)  =  .OOOOOE+00 


CHEM.TIME(YEARS) 


.OOOOOE-kOO 


NO.  MOVES  CX)MPLETED  «  0 


0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 
0  0  0 


00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 


ICONCENTRATION  OF  OXYGEN 

NUMBER  OF  TIME  STEPS  =  0 

TIME(SECONDS)  =  BOOOO 

CHEM.TIME(SECONDS)  =  .OOOOOE+00 
CHEM.TIME(DAYS)  =  .OOOOOB+00 
TIME(YEARS)  =  .OOOOOE+00 
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CHEM.TIME(YEARS) 


X)0000E+00 


NO.  MOVES  COMW.ETED  *  0 


00000000000 

08888888880 

08888888880 

08888888880 

08888888880 

08888888880 

08888888880 

08888888880 

08888888880 

08888888880 

08888888880 

08888888880 

08888888880 

08888888880 

08888888880 

08888888880 

08888888880 

08888888880 

08888888880 

00000000000 


N=  1 

NUMBER  OF  ITERATIONS  =  20 

IHEAD  DISTRIBUTION  -  ROW 

NUMBER  OF  TIME  STEPS  =  1 

TIME(SECONDS)  =  .18935B+09 
TIME(DAYS)  =  .21915E+04 
TIME(YEARS)  =  .60000E+01 

0  .0000000  .0000000  .0000000  .0000000  .0000000  .0000000  BOOOOOO  .0000000  .0000000 

0 
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0  995999998  99.9999998 

99.9999998 

99.9999998 

99.9999998 

99.9999998 

99.9999998  99.9999998 

995999998  0 

0  99.8291477  99.8294635 

993302416 

993310992 

bo 

00 

993310992 

993302416 

995294635  995291477  0 

0  99.6586734  99.6595657 

99.6612187 

993633789 

993649445 

00 

993612187 

99.6595657  995586734  0 

0  99.4866865  99.4880365 

99.4909366 

99.4955451 

993007236 

99.4955451 

99.4909366 

99.4880365  99.4866865  0 

0  99.3131983  99J147908 

993187861 

993270257 

993467482 

993270257 

993187861 

99J147908  993131983  0 

0  99.1386046  99.1399849 

99.1429669 

99.1476276 

99.1528277 

99.1476276 

99.1429669 

99.1399849  99.1386046  0 

0  985622862  98.%32029 

98.9650297 

98.9673115 

98.9689244 

98.9673115 

985650297 

98.9632029  985622862  0 

0  98.7848457  98.7853475 

98.7863989 

98.7874966 

98.7880871 

98.7874966 

98.7863989 

98.7853475  98.7848457  0 

0  98.6071305  98.6074165 

98.6079959 

98.6085203 

98.6087643 

98.6085203 

983079959 

98.6074165  98.6071305  0 

0  98.4289299  98.4290704 

98.4293942 

98.4296452 

98.4297534 

98.4296452 

98.4293942 

98.4290704  98.4289299  0 

0  983503905  98.2504197 

983506138 

983507353 

983507850 

983507353 

983506138 

983504197  983503905  0 

0  98.0720432  98.0720689 

98.0721733 

98.0722224 

98.0722427 

98.0721733 

98.0720689  8.0720432  0 

0  97.8934704  97.8934850 

973932556 

97.8935769 

973935868 

97.8935769 

97.8935556 

97.8934850  973934704  0 

0  97.7145462  97.7144821 

97.7145573 

97.7145736 

97.7145774 

97.7145736 

97.7145573 

97.7144821  97.7145462  0 

0  973361558  97.5361572 

973361729 

973361670 

973361656 

973361670 

973361729 

973361572  973361558  0 

0  973576164  973576574 

97.3576539 

973576438 

97357642 

973576438 

973576539 

97.3576574  973576164  0 

0  97.1784342  97.1783588 

97.1784118 

97.1784250 

97.1784260 

97.1784250 

97.1784118 

97.1783588  97.1784342  0 

0  97.0000002  97.0000002 

97.0000002 

97.0000002 

97.0000002 

97.0000002 

97.0000002 

97.0000002  97.0000002  0 

0  .0000000  .0000000  .0000000  .0000000  .0000000  XXXXXXX) 

.0000000  .0000000  .0000000  0 

IHEAD  DISTRIBUTION  -  ROW 

NUMBER  OF  TIME  STEPS  =  1 

TIME(SECONDS)  =  .18935E+09 

TIME(DAYS)  =  21915E+04 

TIME(YEARS)  =  .60000E+01 

0  00000000000 
0  0  100  100  100  100  100  100  100  100  100  0 

0  0  100  100  100  100  100  100  100  100  100  0 

0  0  100  I  'JC  100  100  100  100  100  100  100  0 

0  0  99  99  99  99  100  99  99  99  99  0 
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0 

0 

99 

99 

99 

99 

99 

99 

99 

99 

99 

0 

0 

0 

99 

99 

99 

99 

99 

99 

99 

99 

99 

0 

0 

0 

99 

99 

99 

99 

99 

99 

99 

99 

99 

0 

0 

0 

99 

99 

99 

99 

99 

99 

99 

99 

99 

0 

0 

0 

99 

99 

99 

99 

99 

99 

99 

99 

99 

0 

0 

0 

98 

98 

98 

98 

98 

98 

98 

98 

98 

0 

0 

0 

98 

98 

98 

98 

98 

98 

98 

98 

98 

0 

0 

0 

98 

98 

98 

98 

98 

98 

98 

98 

98 

0 

0 

0 

98 

98 

98 

98 

98 

98 

98 

98 

98 

0 

0 

0 

98 

98 

98 

98 

98 

98 

98 

98 

98 

0 

0 

0 

98 

98 

98 

98 

98 

98 

98 

98 

98 

0 

0 

0 

97 

97 

97 

97 

97 

97 

97 

97 

97 

0 

0 

0 

97 

97 

97 

97 

97 

97 

97 

97 

97 

0 

0 

0 

97 

97 

97 

97 

97 

97 

97 

97 

97 

0 

0 

0 

0 

0 

0 

0 

0  1 

D  C 

1  0 

0 

0 

1  DRAWDOWN 

00000000000 
00000000000 
0  -99-99-99  -99-99-99-99-99  -99  0 
0-99-99-99  -99-99  -99-99  -99-99  0 
0  -98  -98  -98  -98  -99  -98  -98  -98  -98  0 
0  -98  -98  -98  -98  -98  -98  -98  -98  -98  0 
0  -98  -98  -98  -98  -98  -98  -98  -98  -98  0 
0  -98  -98  -98  -98  -98  -98  -98  -98  -98  0 
0  -98  -98  -98  -98  -98  -98  -98  -98  -98  0 
0  -98  -98  -98  -98  -98  -98  -98  -98  -98  0 
0-97  -97  -97  -97  -97  -97  -97  -97-97  0 
0-97-97  -97  -97  -97-97  -97  -97-97  0 
0-97  -97  -97  -97  -97  -97-97-97  -97  0 
0-97  -97  -97  -97-97  -97-97  -97  -97  0 
0-97  -97  -97  -97-97  -97  -97  -97  -97  0 
0-97  -97  -97  -97-97-97  -97  -97-97  0 
0 -% -% -96 -96 -96 -96 -96 -96 -96  0 
0 -96 -96 -% -96 -% -96 -96 -96 -%  0 
00000000000 
00000000000 

0  CUMULATIVE  MASS  BALANCE  -  (IN  FT**3) 
RECHARGE  AND  INJECTION  =  -37869E+05 
PUMPAGE  AND  E-T  WITHDRAWAL  =  .OOOOOE+00 
CUMULATIVE  NET  PUMPAGE  =  -37869E+05 
WATER  RELEASE  FROM  STORAGE  =  jOOOOOE+00 
LEAKAGE  INTO  AQUIFER  =  .72354E+06 
LEAKAGE  OUT  OF  AQUIFER  =  -.760(r7E+06 
CUMULATIVE  NET  LEAKAGE  =  -36531E+05 

0  MASS  BALANCE  RESIDUAL  =  1338^ 

ERROR  (AS  PERCENT)  =  .17596 
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0  RATE  MASS  BALANCE  -  (IN  Cf^.) 

LEAKAGE  INTO  AQUIFER  »  38213E-02 
LEAKAGE  OUT  OF  AQUIFER  »  •.40142E-02 
NET  LEAKAGE  (QNED  »  -.19293E-03 
RECHARGE  AND  INJECTION  »  .2(X)00E-O3 
PUMPAGE  AND  E-T  WITHDRAWAL  =  jOOOOOE-HX) 
NET  WITHDRAWAL  (TPUM)  *  .20000E-03 


STABILITY  CRITERIA  -  M.O.C. 

FLUID  VELOCITIES 

VMAX^  932E-08  VMAY  =  U6E-06 
VMXBD=  131E-07  VMYBD=  1J29E-06 
TMV  (MAX.  INJ.)  =  .43898E+08 
TTMV  (CELDIS)  =  .19338E+08 
0  TTMV  «  1.93E+07  NTTMV  =  9  NMOV  =  10 
TIM(N)  =  .18935E+09 
TTMEVELO  =  .18935E+08 
TTMEDISP  =  37852E+08 

0  TTMV  =  139E+07  NTTMD  =  2  NMOV  =  lO 


THE  LIMITING  STABILITY  CRITERION  IS  CELDIS 
0  NO.  OF  TIME  STEPS  REQUIRED  TO  COMPLETE  THIS  TIME  STEP  =  20 


ONP 

TIM(N) 

ONPl 

TIM(N) 

ONP 

TIM(N) 

ONPl 

TIM(N) 


=  1709  IMOV  *  1 

=  .18935E+09  TTMV  =  .18935E+08  SUMTCH  *  .18935B+08 
=  1709  IM0V(02)=  1 

=  .18935E+09  TTMV  =  .18935E+08  SUMTCH  =  .18935E+08 
=  1715  IMOV  =  2 

=  .18935E+09  TTMV  =  .18935Bf08  SUMTCH  =  37869E+08 
=  1715  IMOV(02)*  2 

=  .18935E+09  TTMV  =  .18935E+08  SUMTCH  =  37869E+08 
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ONP  »  1715  IMOV  =  3 

TIM(N)  =  .18935B+09  HMV 
ONPl  »  1715  IMOV(02)  = 

TIM(N)  =  .18935E+09  TIMV 
ONP  »  1724  IMOV  =  4 

TTMCN)  =  ,18935E+09  TIMV 
ONPl  =  1724  IMOV(02)  = 

TIM(N)  =  .18935E+09  TIMV 
ONP  =  1740  IMOV  =  5 
TZM(N)  =  J8935B+09  TIMV 
ONPl  =  1740  IMOV(02)  = 

TIM(N)  =  .18935E+09  TIMV 
ONP  =  1748  IMOV  =  6 

TIM(N)  =  J8935E+09  TIMV 
ONPl  =  1748  IMOV(02)  = 

TIM(N)  =  .18935E+09  TIMV 
ONP  =  1757  IMOV  =  7 

TIM(N)  =  .18935E+09  TIMV 
ONPl  =  1757  IMOV(02)  = 

T1M(N)  =  .18935E+09  TIMV 
ONP  =  1757  IMOV  =  8 

TIM(N)  =  .18395E+09  TIMV 
ONPl  =  1757  IMOV(02)  = 

TIM(N)  =  .18935E+09  TIMV 
ONP  =  1760  IMOV  =  9 

TIM(N)  =  .18935E+09  TIMV 
ONPl  =  1760  IMOV(02)=! 

TIM(N)  =  .18935E+09  TIMV 
ONP  =  1764  IMOV  =  10 

TIM(N)  =  .18935E+09  TIMV 
ONPl  =  1764  IMOV(02)  = 

TIM(N)  =  .18935E+09  TIMV 


=>  .18935E-f08  SUMTCH  s  J6804B+08 

3 

s  .18935E408  SUMTCH  >  36804E+08 

=  J8935B+08  SUMTCH  =  .75738E+08 

4 

=  .18935E+OB  SUMTCH  =  .75738B+08 

=  .18935E+08  SUMTCH  =  .94673E+08 

5 

=  .18935E+08  SUMTCH  =  54673E+08 

=  .18935E+08  SUMTCH  =  .11361E+09 

6 

=  .18935E+08  SUMTCH  =  .11 361E+09 

=  .18935E+(»  SUMTCH  =  .13254E+09 

7 

=  .18935E+08  SUMTCH  =  .13254E+09 

=  ,18935E+(»  SUMTCH  =  .15148E+09 

8 

=  .18935E+08  SUMTCH  =  .15148E+09 

=  .18935E+08  SUMTCH  =  .17041E+09 

9 

=  .18935E+08  SUMTCH  =  .17041E+09 

=  .18935E+08  SUMTCH  =  .18935E+09 

10 

=  .18935E+08  SUMTCH  =  .18935E+09 


ICONCENTRATION  OF  CONTAMINANT 

NUMBER  OF  TIME  STEPS  =  1 

DELTA  T  =  .18935E+09 

TIMECSECONDS)  =  .18935E+09 

CHEM.TIME(SECONDS)  =  .18935E+09 
CHEM.TIME(DAYS)  =  21918E+04 
TIME(YEARS)  =  .60000E-t-01 
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CHEM.TIME(YEARS)  »  iSOOOQ&fOl 


NO.  MOVES  COMPLETED  =  10 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

30 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

35 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

30 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

26 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

18 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

9 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

CHEMICAL  MASS  BALANCE 

MASS  IN  BOUNDARIES  =  .0000(£-K)0 
MASS  OUT  BOUNDARIES  =  -.82523E-03 
MASS  PUMPED  IN  =  37869E+07 
MASS  PUMPED  OUT  =  DOOOOE+OO 
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MASS  LOST  W.  BIODEO.  »  .7897411+06 
MASS  LOST  BY  RADIO.  DCY«  BOOOOB+OO 
MASS  LOST  BY  ANAER.  DCYs  .0000(£+00 
MASS  LOST  BY  REAER.  DCY»  XXMOOE+OO 
MASS  ADSORBED  ON  SOLIDS^  DOOOIffi+OO 

INITIAL  MASS  ADSORBED  »  DOOOOE+OO 
INFLOW  MINUS  OUTFLOW  =  37869E+07 
INITIAL  MASS  DISSOLVED  »  OOOOOE+00 
PRESENT  MASS  DISSOLVED  =  32136E+07 
CHANGE  MASS  DISSOLVED  =  .40033E+07 
CHANGE  TOTLA4ASS  STORED*  .40033E+07 


COMPARE  RESIDUAL  WITH  NET  FLUX  AND  MASS  ACCUMULATION: 

MASS  BALANCE  RESIDUAL  =  -31638E+06 
ERROR  (AS  PERCENT)  =  -37140E+OI 


ICONCENTRATION  OF  OXYGEN 

NUMBER  OF  TIME  STEPS  =  1 

DELTA  T  =  .18935E+09 

TIME  (SECONDS)  =  .18935E+09 

CHEM.  TIME  (SECONDS)*  .18935E+09 
CHEM.  TIME  (DAYS)*  31915E+04 
TIME  (YEARS)  =  .60000E+01 

CHEM.  TIME  (YEARS)*  .60000E+01 
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NO.  MOVES  COMPLETED  =  10 

00000000000 

08888888880 

08888888880 

08888788880 

08888088880 

08887078880 

08882028880 

08880008880 

08880008880 

08881018880 

0888S0S8880 

08887078880 

08888688880 

08888888880 

08888888880 

08888888880 

08888888880 

08888888880 

08888888880 

00000000000 

CHEMICAL  MASS  BALANCE  FOR  OXYGEN 

MASS  IN  BOUNDARIES  =  .57883E+07 
MASS  OUT  BOUNDARIES  *  -.60806E+07 
MASS  PUMPED  IN  =  .OOOOOE+00 
MASS  PUMPED  OUT  =  .OOOOOE+00 
MASS  LOST  W.  BIODEG.  =  23692E+07 
INFLOW  MINUS  OUTFLOW  =  -.20225E+06 
INITIAL  MASS  DESOLVED  =  .24300E+08 
PRESENT  MASS  DESOLVED  =  21722E+08 
CHANGE  MASS  DISSOLVED  =  -.20918E406 
CHANGE  TOTL.MASS  STORED=  -20918E+06 


COMPARE  RESIDUAL  WITH  NET  FLUX  AND  MASS  ACCUMULATION 
FOR  OXYGEN: 

MASS  BALANCE  RESIDUAL  =  -.83069E+05 
ERROR  (AS  PERCENT)  =  -.14351E+01 

COMPARE  INITIAL  MASS  STORED  WITH  CHANGE  IN  MASS  STORED 
FOR  OXYGEN: 

ERROR  (ASPERCWD  =  33778E+00 
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0  TIME  VERSUS  HEAD  AND  CONCENTRATION  AT  SELECTED 
OBSERVATION  POINTS 

PUMPING  PERIOD  NO.  1 


0  STEADY-STATE  SOLUTION 


OBS.WELL  NO. 


X 


Y 


N 


HEAD  (FT)  CONC  (mg/1)  TIME  (Yrs) 


1  6  10 


0 

98.6 

.0 

.000 

1 

98.6 

.0 

.600 

2 

98.6 

.0 

UOO 

3 

98.6 

.0 

1300 

4 

98.6 

.0 

2.400 

5 

98.6 

.0 

3.000 

6 

98:6 

.0 

3.600 

7 

98.6 

2:6 

4200 

8 

98.6 

9.7 

4300 

9 

98.6 

13.6 

5.400 

10 

98.6 

173 

6.000 

1  Bioplume  II  -  Sample  Data 


0  TIME  VERSUS  HEAD  AND  C0NCENTRATI0N(02)  AT  SELECTED 
OBSERVATION  POINTS 

PUMPING  PERIOD  NO.  1 


STEADY-STATE  SOLUTION  FOR  OXYGEN 
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OBS.WELL  NO. 

1 


Y 

N 

HEAD  (FT) 

CONC  (mg/l) 

10 

1 

0 

98.6 

jO  8D 
8j0 

iSOO 

2 

98.6 

8j0 

1200 

3 

98.6 

8j0 

1.800 

4 

98.6 

7.7 

2.400 

5 

98.6 

5.6 

3.000 

6 

98.6 

5 

3.600 

7 

98.6 

.0 

4.200 

8 

98.6 

.0 

4.800 

9 

98.6 

jO 

5.400 

10 

98.6 

.0 

6IX)0 

TIME(Yrs) 


jOOO 
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Appendix  D 
Data  Input  Fbnnats 


Appendix  D:  Data  Input  Formats 


Card 

Colunm 

Format 

Variable 

Definition 

1 

1-80 

10A8 

TITLE 

Description  of  problem 

2 

1-4 

14 

NTTM 

Maximum  number  of  time 
steps  in  a  pumping  period  (1..100)*. 

5-8 

14 

NPMP 

Number  of  pumping  periods.  If 
NPMP  >  1,  then  data  set  11  must 
be  completed. 

100 


9-12 

14 

NX 

Number  of  nodes  in  the  x 
direction  (3.  JO). 

13-16 

14 

NY 

Number  of  nodes  in  the  y 
directitm  (3.30). 

17-20 

14 

NPMAX 

Maximum  number  of  particles 
(1.3100). 

21-24 

14 

NPNT 

Time  step  interval  pinting 
hydraulic  and  ctonical  output  data 
(1-). 

25-28 

14 

NTTP 

Number  of  iteration  parametmx 
(usually  4..7). 

29-32 

14 

NUMOBS 

Number  of  observation  points  to 
be  specified  in  a  following  date 
set  (0.3). 

33-36 

14 

ITMAX 

Maximum  allowable  number  of 
iterations  in  ADD*  (usually  100.  JOO). 

37-40 

14 

NREC 

Number  of  pumping  or  injection 
wells  to  be  ^tecified  in  a 
following  date  set  (0.30). 

*  The  possible  ranges  which  a  value  may  take  are  indicated  in  parentheses.  Theform  (X..Y)  indicates 
a  lower  bound  and  Y  is  an  upper  bound.  The  form  (X.«)  indicates  that  X  is  a  lower  bound  and  there  is 
bound. 

Card  Column 

Format 

Variable 

Definition 

2  41-44 

14 

NPTPND 

Initial  number  of  particles  per 
node  (433.9). 

45-48 

14 

NCODES 

Number  of  node  identification 
codes  to  be  ^tecified  in  a 
following  date  set  (1..9). 

49-52 

14 

NPNTMV 

Particle  movement  interval 
(IMOV)  for  printing  chemical 

that  X  is 
no  upper 
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output  data  (0  to  print  tmly  at 
end  of  time  steps). 


53-56 

14 

NPNTVL 

Option  for  printing  computed 
velocities  (Ossdo  not  print; 
l=print  for  first  time  step; 
2sprint  for  all  time  steps). 

57-60 

14 

NPNTD 

Option  for  printing  computed 
di^rsion  equation  coefficients 
(0=do  not  print;  l=print  for  first 
time  step;  2=print  for  all  time 
steps). 

61-64 

14 

NPDELC 

Option  for  printing  computed 
changes  in  concentration 
(0=do  not  print;  1  =print). 

65-68 

14 

NPNCHV 

Option  to  punch  velocity  data 
(Osdo  not  punch;  1  =punch  the 
node  velocities  on  unit  7). 

69-72 

14 

NREACT 

Option  for  retardation  and 
decay  (0=retardation  factor  is 
equal  toland  no  decay; 
l=retardation  factor  > land/or 
decay). 

Card  Column 

Format 

Variable 

Definition 

3  1-5 

05.0 

PINT 

Pumping  period  in  years 
(0.01. .99.99). 

6-10 

G6.0 

TOL 

Convergence  criteria  in  ADD* 
(usually  <  0.01). 

11-15 

G5.0 

POROS 

Effective  porosity  (0.01. .1). 
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Card 


4 


16-20 

G5.0 

BETA 

Characteristic  Imgth 
(Iragitudinal  dispersivity)  in 
feet  (0X)1. 59.99). 

21-25 

G5.0 

S 

Storage  coefficient  (set  8=0  for 
steady  flow  problems). 

26-30 

G5.0 

TEMX 

Tune  increment  multiplier  for 
transient  flow  problems 
(0.01.59.99).  TTMX  is 
disregarded  if  S=0. 

31-35 

G5.0 

TINIT 

Size  of  initial  time  step  in 
seconds  (0...).  UNIT  is 
disregarded  if  S=0. 

36-40 

G5.0 

XDEL 

Width  of  finite  difference  cell 
in  the  x  direction  in  feet 
(0.1  ..999.9). 

41-45 

G5.0 

YDEL 

Width  of  finite  difference  cell  in 
the  y  direction  in  feet 
(0.1  ..9995). 

46-50 

G5.0 

DLTRAT 

Ratio  of  nansverse  to 
longitudinal  dispersivity 
(0.001  ..1). 

51-55 

G5.0 

CELDIS 

Maximum  cell  distance  per 
particle  move  (0.001. .1). 

56-60 

G5.0 

ANFCTR 

Ratio  of  Tyy  toTxx 
(0.001.5.99). 

Colunm 

Format 

Variable 

Definition 

Free  Format 

DK 

Distribution  coefficient  (L^  /M). 

(separated 
by  spaces 

RHOB 

Bulk  density  of  the  solid  (M/L^ 

or  commas) 

THALF 

Half-life  of  the  solute  (seconds) 

DECl 

Anaerobic  decay  coefficient 
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DEC2 


Reaenticm  coefficient  (day'* ) 


Data  Number  Format  Variable  Definition 

Set  of  Lines 


1  NUMOBS  212  DCOBS,  Coordinates  of  observatitm 

lYOBS  points.  This  data  set  is  not 

used  if  NIJMOBS=0. 


2  NREC  212^08.2 


DC.  lY 


REC 

CNRECH 


Coordinates  of  pumping  (■*■)  or 
injection  (-}  wells  for 
contaminant  or  oxygenated 
water 

Pumping/injection  rate  in  cfs 

Concentration  of  injected 
contaminated  water 


CNRECO  Concentration  of  injected 

oxygenated  water 


Data  Number  Format  Variable  DeHnition 

Set  of  Lines 


3  1  or  NY  II  INPUT  Parameter  card  for 

transmissivity 

(Osconstant  transmissivity  is 
defined  by  FCTR; 
l=transmissivity 
is  read  from  following  array). 


G10.0 

FCTR 

Constant  transmissivity  in  ft  /s 
OR  factor  to  multiply 
transmissivity  array. 

20G4.  1 

VPRM 

Array  for  temporary  storage  of 
transmissivity  data  in  ft*  /s.  For 

an  anisotropic  array,  enter  the 
values  of  Txx  and  the  values 

for  Tyy  will  be  computed  by 
multiplying  by  ANFCTR. 

4  1  or  NY  II  INPUT 


G10.0  FCTR 


20G3.0  THCK 


Parameter  card  for  thickness 
(0=constant  thickness  is 
defined  by  FCTR;  l=thickness 
is  read  from  following  array). 

Constant  thickness  in  feet,  OR 
factor  to  multiply  thickness 
array. 

Array  of  saturated  thickness  in 
feet. 


Data  Number  Format  Variable  Defmition 

Set  of  Lines 


5  1  or  NY  II  INPUT 


G10.0  FCTR 


Parameter  card  for  recharge 
(0=constant  recharge  is 
defined  by  FCTR;  l=recharge 
is  read  from  following  array) 

Constant  diffuse  recharge  (-) 
or  discharge  (-f)  in  ft/s  OR 
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factor  to  multiply  recharge 
array 

20G4.1  RECH  Array  of  diffuse  recharge  (-)  or 

discharge  (+)  in  ft/s. 


6  1  or  NY  II  INPUT  Parameter  card  for  node 

identification  (Osall  nodes 
identified  by  PCTR;  I-node 
identificatimis  in  following 
array) 

G10.0  FCTR  Node  identification  OR  factor 

to  multiply  node  identification 
array 

2011  NODEID  Node  identification  matrix 

(used  to  define  constant-head 
nodes  or  other  boundary 
conditions  and  stresses). 


Data  Niunber  Format  Variable  Definition 

Set  of  Lines 

7  NCODES  12  ICODE  Instructions  for  using  the 

NODEID  array.  When  NODE 
ID=  ICODE,  then  the  following 
factors  are  set.  Otherwise,  the 
values  remain  set  as  they  were 
previously. 

4G  102  FCTRl  Leakance 
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FCTR2 

Crmcentratitm  of  contaminated 
water 

FCTR4 

Concentration  of  oxygenated 
water 

FCTR3 

Diffuse  recharge  (-)  or 
discharge  (-•■) 

12 

OVERRD 

If  OVERRD=0,  then  the  value 
of  RECH  is  not  changed.  If 
OVERRD  is  nonzero,  then  the 
value  of  RECH  is  set  to  FCTR3 

8  lor  NY  11  INPUT 


G10.0  FCTR 


2004.0  WT 


Parameter  card  for  water  table 
(0=constant  water  table 
defmed  by  FCTR;  l=water 
table  is  read  from  following 
array) 

Initial  water  table,  potentio- 
metric  elevation,  or  constant 
head  in  stream  or  source  bed 
in  feet  OR  factor  to  multiply 
water  table  array. 

Array  of  initial  water  table, 
potentiometric  elevation,  or 
constant  bead  in  stream  or 
source  bed  in  feet. 


Data 

Number 

Format 

Variable 

Definition 

Set 

of  Lines 

9  1  or  NY  II  INPUT  Parameter  card  for  initial 

contaminant  concentration 
(0=constant  concentration 
defmed  by  FCTR; 
l=contaminant  concentration  is 
read  from  following  array) 
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G10.0  FCTR 


20G4i)  CONC 


Initial  contaminant 
concentratitm  in  aquifer  OR 
factor  to  multiply  contaminant 
concentration  array 

Array  of  initial  contaminant 
concentration  in  aquifer. 


10  1  or  NY  II  INPUT 


G10.0  FCTR 


20G4.0  CONC  1 


Parameter  card  for  initial 
oxygen  concentration 
(D=constant  concentration 
defined  by  FCTR;  l=oxygen 
concentration  is  read  from 
following  array) 

Initial  oxygen  concentration  in 
aquifer  OR  factor  to  multiply 
oxygen  concentration  array 

Array  of  initial  oxygen 
concentration  in  aquifer. 


Data  Number  Format  Variable  Definition 

Set  of  Lines 


11  1  II  ICHK  Parameter  to  check  whether 

any  revisions  are  desired 
(l=revision  is  desired,  more 
data  to  follow;  0=no  revision 
desired,  end  of  data  set).  This 
data  set  allows  13  timing, 
printing,  and  pumping 
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NREC 


1014. 

NTIM 

variables  to  be  revised  fm 
each  pumping  period.  Data  set 
1 1  can  only  be  used  if 

NPMP  >  1.  The  sequence  of 
cards  in  this  data  set  must  be 
repeated  NPMP-1  times  (for 
each  pumping  period  after  the 
first). 

Previously  defmed  variables 

3G5.0 

NPNT, 

for  cards  2  and  3  which  will 

NITP. 

be  revised  for  the  next 

ITMAX, 

pumping  period.  This  card  is 

NREC. 

used  only  if  ICHK=1. 

212, 

NPNTMV. 

NPNTVL. 

NPNTD, 

NPDELC. 

NPNCHV. 

P  I  NT. 

TI  MX. 

T  I  N  IT 

IX.  lY,  REC. 

Previously  defined  variables 

3G8.2 

CNRECH. 

for  data  set  2  which  will 

CNRECO 

be  revised  for  the  next 

pumping  period.  This  card  is 
used  only  if  ICHK*1  and 
NREC>0 
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Appendix  E;  BIOn.UME  11™  Sensitivity  Analysis 


Sensitivity  Analysis 

In  order  to  define  which  parameters  have  the  mcMt  effect  m 
Inodegradation  in  BIOn^UME  U™,  the  following  detailed  sensitivity  analysis  was 
perftnnned.  The  parameters  that  were  investigated  included:  hydraulic 
conductivity,  diqpersivity  (Icmgitudinal  and  transverse),  porosity,  reaeratitm,  and 
retardation. 

A  hypothetical  ctmtaminant  plume  was  generated  using  a  single  ctmtinuous  hydrocarbtm  source.  The  above 
mentioned  parameters  were  then  varied  individually  to  tktermine  their  effect  on  biodegradatitm.  The  results  from 
the  sensitivity  analyses  indicate  that  biodegradation  in  the  model  is  most  sensitive  to  hydraulic  omductivity,  the 
coefficient  of  reaeration,  and  the  coefficient  of  anaerolnc  decay.  The  following  input  data  was  utilized  in  the 
base  run: 


Simulation  time  10  years 

Grid  size  20  x  30 

Cell  size  SO  ft  x  SO  ft 

Porosity  OJ 

Longitudinal  Dispersivity  10  ft 

Transverse  Dispersivity  3  ft 

Txx  .002S  ft  /s 

Tyy  .002S  ft  /s 

Aquifer  thickness  25  ft 

Hydraulic  Gradient  4J9E-3  ft/ft 

Injection  well  at  cell  X=10,Y=10 

Injection  rate  0.0002  cfs 

Cone,  of  contaminant  in  injected  water  ISO  mg/1 

Cone,  of  oxygen  in  injected  water  0.0  mg/1 

Initial  Cone,  of  oxygen  8D  mg/1 

Cone,  of  natural  recharge  of  oxygen  8X)  mg/I 


A  detailed  disctission  of  the  sensitivity  analyses  is  included  in  the  following 
sections. 

23.1  Variation  of  Concentrations  with  Hydraulic  Conductivity 

The  hydraulic  conductivity  (K)  was  varied  from  10^  ft/sec  to  10*^  fr/sec. 

Hgure  23  is  a  plot  of  the  contaminant  and  oxygen  concentrations  along  the  centerline  of  the  plume  (cross 
section  A  -  A,  Rgure  22)  for  three  values  of 

hydraulic  conductivity.  It  can  be  seen  that  the  hydraulic  conductivity  has  a 
significant  effect  on  biodegradation.  The  maximum  contaminant  concentration 

varied  from  29.6  mg/1  (K  =  10^  ft/sec)  to  130  mg/1  (K  =  fr/sec).  The  change  in  biodegraded  mass  with 

hydraulic  conductivity  is  illustrated  in  Figure 

2.4. 


Ill 


2  J.2  VariatioB  of  CoBcentratloBS  with  ReUrdatloB 


The  effect  of  retardation  on  biodegradation  was  studied  by  using  a 
retardation  factor.  R,  greater  than  1.  Figure  15  presents  the  variation  in 
contaminant  and  oxygen  concentrations  along  the  centerline  of  the  plume 
(sectitm  A  -  A,  Fgure  12)  for  three  values  of  R.  It  can  be  seen  that  the  mass 
of  hydrocarbon  remaining  at  the  end  of  the  simulation  period  decreases  with 
increasing  values  of  R.  The  percent  of  mass  biodegraded  relative  to  the  total 
dissolved  mass,  however,  decreases  for  increasing  values  of  R  (Hgure  2.4b). 

Table  2.1  lists  the  percent  of  mass  biodegraded  relative  to  the  total 
dissolved  mass  and  the  percent  of  mass  biodegraded  relative  to  the  total 
stored  mass.  Table  2.1  also  lists  the  percent  of  mass  adsorbed  relative  to  the 
total  stored  mass  for  the  three  values  of  R.  Table  2.1  indicates  that  the 
percent  of  biodegraded  mass  decreases  with  increasing  values  of  retardation 
(columns  A  and  B). 

Table  2.1  -  Percent  Biodegraded  Mass  as  a  Function  of 


the  Retardation  Factor  (R) 


R 

A 

B  C 

1 

30.25 

30.25  0.00 

2 

27.73 

13.86  50D0 

3 

25.74 

838  6637 

A  =  BM/TDM 

B  *  BM/TSM 

C  =  AM/TSM 

BM  =  Biodegraded  Mass 

AM  =  Adsorbed  Mass 

TDM  s  Dissolved  Mass  -4-  Biodegraded  Mass 
TSM  s  Dissolved  Mass  +  Biodegraded  Mass 
^Adsorbed  Mass 

23.3  Variation  of  Concentrations  with  Dispcrsivity 

The  variation  of  contaminant  concentrations  with  dispersivity  was 
examined  by  looking  at  the  longitudinal  and  transverse  dispersivities 
independently. 


233.1  Longitudinal  Dispcrsivity 

The  longitudinal  dispersivity  was  varied  from  10  ft  to  1  ft.  Fgure  2.6 
[H'esents  the  variation  in  contaminant  and  oxygen  concentrations  along  the 
centerline  of  the  contaminant  plume  (section  A  -  A,  Figure  2.2)  for  three 
values  of  longitudinal  dispersivity.  It  can  be  seen  that  the  longitudinal 
dispersivity  also  has  a  slight  effect  on  biodegradation.  The  maximum 
contaminant  concentrations  varied  from  29.6  mg/1  (10  ft)  to  41.7  mg/1  (1  ft). 
The  change  in  biodegraded  mass  with  longitudinal  dispersivity  is  illustrated  in 
Figure  2.7. 
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2  J  J.2  TraasverM  Dispersivity 

The  transverse  diqwrsivity  was  varied  from  1  ft  to  5  ft.  Figure  2.8  shows 
the  variation  of  contaminant  and  oxygen  concentrations  with  transverse 
dispersivity  along  a  transverse  cross  section  through  the  centerline  of  the 
plume  (cross  section  B  -  B,  Figure  2.2).  The  transverse  dispersivity  does  not 
seem  to  have  an  appreciable  effect  on  biodegradation.  The  areal  extent  of  the 
plume  is  not  very  sensitive  to  the  transverse  dispersivity,  however,  the 
maximum  concentrations  exhibit  a  wide  range  of  variation.  The  maximum 
contaminant  concentrations  varied  from  26.3  mg/1  (  5  ft)  to  3S.S  mg/1  (1  ft). 

The  change  in  biodegraded  mass  with  transverse  dispersivity  is  illustrated  in 
Figure  2.7 

2  J.4  Variation  of  Concentrations  with  PorMity 

The  porosity  was  varied  from  0.23  to  0.7.  Figure  2.9  shows  the  variation 
of  contaminant  and  oxygen  concentrations  with  porosity  along  the  centerline 
of  the  plume  (section  A  -  A,  Figure  2.2).  It  is  evident  that  porosity  does  not 
have  a  significant  effect  on  biodegradation.  The  maximum  contaminant 
concentrations  varied  from  29.0  mg/1  (  n  =  0.23)  to  33.0  mg/1  (n  =  05).  The 
change  in  biodegraded  mass  with  porosity  is  illustrated  in  Fgure  2.4b. 

23.5  Variation  of  Concentrations  with  Rcaeratlon  Coefficient 

The  reaeration  coefficient,  k,  was  varied  from  0.0  day  to  0.003  day  . 

Figure  2.10  presents  the  variation  of  contaminant  and  oxygen  concentrations 
with  k  along  the  centerline  of  the  plume  (cross  section  A  -  A,  Fgure  22).  It  is 
evident  that  the  coefficient  of  reaeration  has  a  significant  effect  on 
biodegradation.  The  areal  extent  of  the  contaminant  plume  as  well  as  the 
maximum  concentrations  exhibit  a  wide  range  of  variation  with  k.  The 
maximum  concentrations  varied  from  29.6  mg/1  (k  =  0.0  day  )  to  17.6  mg/I  (k  = 
0.003  day  ).  The  change  in  biodegraded  mass  with  k  is  presented  in  Hgure 
2.4a. 

2.4  Model  Output 

Typical  output  from  BIOPLUME  11  includes  an  oxygen  and  hydrocarbon 
distribution  matrix  at  selected  points  in  time.  These  matrices 
can  be  plotted  as  contour  plots  (Hgure  4.1)  or  surface  plots  using  the  SURFER 
graphics  package  (Golden  Software,  1987).  The  grafdiics  option  in  the 
BIOPLUME  II  preprocessor  will  transform  the  oxygen  and  hydrocarbon  matrices 
to  the  required  format  for  direct  use  in  SURFER. 

It  can  be  seen  from  Figure  4.1  that  whenever  hydrocarbon  is  present  in 
relatively  high  concentrations,  then  oxygen  is  absent.  The  oxygen  plume  forms 
an  envelope  for  the  hydrocarbon  plume  with  oxygen  concentrations  gradually 
increasing  to  initial  background  levels  as  one  moves  away  from  the 
contaminant  plume.  The  model  output  also  includes  a  mass  balance 
computation  for  oxygen  and  hydrocarbon  at  the  selected  points  in  time.  The 
dissolved  mass  present  in  the  system  for  each  is  computed,  as  well  as  the 
biodegraded  mass.  The  hydrocarbon  mass  balance  computation  details  the 
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biodegraded  mass  due  to  the  different  processes  availaUe  in  the  model 
(aerobic,  anaerobic,  reaeration  and  radioactive  decay).  It  is  noted  at  this  point 
that  a  detailed  analysis  of  mass  balance  errors  computed  in  BIOPLUME  1 1  is 
being  performed  for  a  variety  of  conditimis  and  geometries.  The  results  of  the 
analysis  will  be  included  in  future  updates  to  the  manual. 

Samasary 

BIOPLUME  11  simulates  hydrocarbon  transport  under  the  influence  of 
oxygen  limited  biodegradation.  A  dual  particle  mover  concept  is  used  to 
cmnpute  an  oxygen  plume  and  a  hydrocarbon  plume.  An  instantaneous  reaction 
between  the  solute  (hydrocarbon)  and  the  sub^te  (oxygen)  is  assumed  and 
the  method  of  superposition  is  utilized  to  represent  the  reaction  between  the 
two.  An  independent  mass  balance  is  performed  for  oxygen  and  hydrocarbon 
and  is  adjusted  to  account  for  the  mass  loss  due  to  biodegradation. 

The  model  can  be  used  to  simulate  naturally  occuring  biodegradation 
processes  and  to  simulate  in-situ  restoration  processes.  Injection  wells  can  be 
used  as  oxygen  sources  in  the  model.  Three  other  sources  of  oxygen  are 
included  in  BIOPLUME  1 1:  (1)  dissolved  oxygen  in  the  aquifer;  (2)  natural 
recharge,  and;  (3)  oxygen  exchange  fron  the  unsaturated  zone. 

The  biodegraded  mass  in  the  model  is  most  sensitive  to  hydraulic 
conductivity,  the  coefficient  of  reaeration,  and  the  coefficient  of  anaerobic 
decay.  The  model  has  been  applied  to  two  sites:  a  wood  creosoting  process 
waste  site  in  Conroe.  Texas  (Borden  et  al.,  1986)  and  a  jet  fuel  spill  site  in 
Traverse  City,  Michigan.  The  model  application  to  the  Traverse  City  site  is 
presently  being  submitted  for  publication.  The  model  provided  a  good  nuitch 
to  field  conditions  at  both  sites.  BIOPLUME  11  is  presently  being  used  to 
design  an  in-situ  bioreclamation  field  experiment  at  the  Traverse  City  field 
site.  Results  from  the  experiment  will  also  be  published  in  the  literature. 
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FIGURE  2.1  -  Simplified  Flowchart  of  BIOPLUME  11 


x>axi> 


y-axis 

Section  A-A  is  the  Centerline  of  the  Plume 


Section  B*6  is  the  Trsnsverse  Section  of  the  Plume 


FIGURE  2.2  •  Schematic  of  the  Centerline  and 
Transverse  Section  of  a  Plume 
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Figure  2.3  —  Concentration  Distributions  for 
Various  Values  of  Hydraulic  Conductivity 


117 


Figure  2.4a  —  Variation  of  Biodegraded  Mass 

with  Various  Parameters 
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Rgure  2.4b  —  Variation  of  Biodegraded  Mass 

with  Various  Parameters 
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Oxygen  (mg/I)  Contaminant  (mg/I) 


Distance  Along  Centerline  of  Plume  (feet) 


Distance  Along  Centerline  of  Plume  (feet) 


Rgure  2.5  —  Concentration  Distributions  for 
Various  Values  of  Retardation 
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Figure  2.6  —  Concentration  Distributions  for 
Various  Values  of  Longitudinal  Dispersivity 
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Trans.  Dispersivity,  /?*  (ft)  Long.  Dispersivity.  fi  (ft) 


Figure  2.7  —  Variation  of  Biodegraded  Mass 

with  Dispersivity 
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Oxygen  (mg/I)  Contaminant  (mg/I) 


Distance  Along  Transverse  Section  of  Plume  (ft) 


Distance  Along  Transverse  Section  of  Plume  (ft) 


Figure  2.8  —  Concentration  Distributions  for 
Various  Values  of  Transverse  Dispersivity 
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Oxygen  (mg/I)  Contaminant  (mg/I) 


Distance  Along  Centerline  of  Plume  (ft) 


Distance  Along  Centerline  of  Plume  (ft) 


Figure  2.9  —  Concentration  Distributions  for 
Various  Values  of  Porosity 
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Oxygen  (mg/I)  Contaminant  (mg/I) 


Distance  Along  Centerline  of  Plunne  (feet) 


Figure  2.10  —  Concentration  Distributions  for 
Various  Values  of  the  Reaeration  Coefficient 
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